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By Franklin W. Diederich and Martin Zlotnick 


SUMMARY 


Spanwise lift distributions have been calculated for nineteen 
unswept wings with various aspect ratios and taper ratios and with a 
variety of angle -of -attack or twist distributions, including flap and 
aileron deflections, by means of the Weissinger method with eight control 
points on the semispan. Also calculated were aerodynamic influence coef- 
ficients which pertain to a certain definite set of stations along the 
span, and several methods are presented for calculating aerodynamic 
influence coefficients for stations other than those stipulated. 

The information presented herein can be used in the analysis of 
untwisted wings or wings with known twist distributions, as well as in 
aeroelastic calculations involving initially unknown twist distributions. 


INTRODUCTION 


In the design and development of an airplane, a knowledge of the 
spanwise lift distribution on the wing is important in predicting the 
structural loads and the stability characteristics. For high-speed air- 
planes having flexible wings, the calculation of the spanwise lift dis- 
tribution is an aeroelastic rather than a purely aerodynamic problem. 

In aeroelastic calculations means are required for calculating the span- 
wise lift distribution for angle -of -attack (or twist) distributions which 
are initially unknown. Aerodynamic influence coefficients constitute the 
most convenient of these means. 

One of the most satisfactory techniques developed in recent years 
for calculating the spanwise lift distribution on a wing in subsonic flow 
has been the Weissinger L-method (ref. l), which can be applied to a 
large variety of plan forms and yields solutions of sufficient accuracy 
for all practical purposes without requiring an unduly long time for the 
calculations. This method may be considered as a simplified lifting- 
surface theory because the calculation of the lift on the wing is treated 
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as a boundary -value problem, the boundary condition being that the down- 
wash angle induced by the bound and trailing vortices is equal to the 
geometric angle of attack at the three -quarter -chord line. 

In the present paper, symmetrical and antisymmetrical lift distri- 
butions and some associated aerodynamic parameters have been calculated 
by means of the Weissinger method with eigjit control points on the semi- 
span for several continuous and discontinuous angle -of -attack conditions 
on nineteen unswept wings having various aspect ratios and taper ratios. 

A convenient matrix formulation of the Weissinger method was used in con- 
junction with the Bell Telephone Laboratories X-66744 relay computer at 
the Langley Laboratory to make the calculations. This formulation is 
described in appendix A. 

Aerodynamic influence coefficients have been calculated for these 
nineteen wings for a certain prescribed set of stations along the span 
and are presented herein, and several methods for calculating aero- 
dynamic influence coefficients for any arbitrary set of stations from 
the numerical results of this paper are also presented. The influence 
coefficients calculated in this manner can be used in aeroelastic anal- 
yses similar to that of reference 2. 


SYMBOIS 


A 

b 


4>ail 


b f 

C BM 


% 

Cl 1/2 


aspect ratio 
wing span 
aileron span 
flap span 

root bending-moment coefficient for unit angle of attack, 
2 x Bending moment 

qSb 

induced-drag coefficient at a unit angle of attack 


lift coefficient at a unit angle of attack 
lift-curve slope per radian for additional-type loading 


lift coefficient for half of antisymmetrically loaded wing at 


a unit tip angle of attack, 


: L l/2 
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C 

C 

C 

C 


rolling -moment coefficient, 


Rolling moment 
qSb 


Id = " C * 


P 



coefficient of damping in roll 

rolling -moment coefficient for -unit aileron deflection 


c 

c 

c l 

|_ Ic bmJ 

L Ic lJ 

LVI 

L Ic iJ 

L l/2 

M 

M 

q 

s 

v 


y 

a 


wing chord 

average chord, s/b 

section lift coefficient 

integrating matrix for 
integrating matrix for 
integrating matrix for 
integrating matrix for 
lift on one semispan 


Crj 4 (see appendix A) 


C-^ (see appendix A) 


C^l/2 ^ see a PP® n< *i x 

Ci (see appendix A) 


A) 


free-stream Mach number 
aerodynamic-influence -coefficient matrix 
dynamic pressure 
wing area 

free-stream velocity 
lateral ordinates 

spanwise location of discontinuity 


spanwise center-of -pressure location 
angle of attack, radians 




k 

— 
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% 

dcj/d& 

effective angle of attack for unit flap deflection, r— 

dc ^ /dQ£ 


r 

vortex strength 



r* 

kr * 

dimensionless vortex strength, — = c cj 

bV 



5 

flap or aileron deflection angle, radians 



e = 

cos -1 y* 



0 O = 

-1 * 
cos A y Q 



■a s 

-1 * 
COS T| 



A' 

sweepback angle, deg 



X 

taper ratio 



Subscripts: 



a 

ant i symme t r i c al 



ail 

aileron 



C 

continuous 



D 

discontinuous 



f 

flap 



L 

left 



R 

right 



s 

symmetrical 



t 

tip 



Superscript: 



* 

dimensionless with respect to semispan b/2 
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Matrix notation: 

[_ | row matrix 

column matrix 


{} 

[] 

LI 

L 1 ! 


general matrix (not a row or a column matrix, but need not 
be square ) 

diagonal matrix 

unit (identity) matrix 

In matrix notation, a prime indicates the transpose of the matrix. 

PRESENTATION OF CALCULATED RESULTS 
Spanwise Lift Distributions 


Geometric characteristics of the nineteen plan forms treated in 
this paper are indicated in table I. Lift distributions due to the 
following continuous symmetric and antisymmetric angle -of -attack dis- 
tributions have been calculated for each of these plan forms: 

Symmetric angle -of -attack distributions: 

Constant (a = l) 

Linear (a = | y* | ) 

Quadratic (a = y*^) 

Cubic (a = | y*^ | ) 

Straight-line (a = -^r|y*|) 

Antisymmetric angle -of -attack distributions: 

Linear (a = y* ) 

Quadratic (a = y*^ for y* ^ 0; a = -y*^ for y* ^ o) 

Cubic (a = y*^) 

Quartic (a = y*^ for y* ^ 0; a = -y*^ for y* % o) 

Qulntic (a = y* 5 ) 
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The straight-line angle -of -attack condition was included to repre- 
sent actual structural twists where the surface of the wing is generated 
by straight lines so that the product c*oc, the deflection of the leading 
edge, varies linearly with y*; that is, 

c*a = c t *y*a t 


or, for unit twist at the tip, 



For untapered wings, the straight-line lift distribution is the same as 
the linear lift distribution, and for wings of zero taper ratio, it is 
undefined. 


Lift distributions for flap -type and aileron -type angle -of -attack 
distributions are also presented. A correction which has been made for 
the spanwise discontinuity in the angle of attack is derived in 
appendix B. The values of b f /b and b ail /b (ratios of the flap span 

to the total span and the aileron span to the total span, respectively ) 
for which the lift distributions have been calculated are 0.1, 0.2, 0.3> 
0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0. As is usual, the flaps have been 
taken to be inboard and the ailerons outboard. The lift distribution 
for any flap or aileron configuration may be obtained, however, by 
linear superposition; thus, the lift distribution for an outboard flap 
extending, for example, from y* = 0.5 to y* =1.0 can be .obtained 

by subtracting the lift distribution for the inboard flap = 0 . 5 ) 


from the additional lift distribution 



A similar procedure 


can be used for inboard ailerons. 


The lift distributions pertaining to each plan form are given in 
one figure, parts (a) and (b) showing the lift distribution due to sym- 
metric and antisymmetric continuous angle -of -attack distributions, 
respectively, and parts (c) and (d) showing the lift distribution due 
to flaps and ailerons, respectively. Lift distributions for plan forms 
with an aspect ratio approaching zero have been taken from reference 3 
and are Included herein in figure 1 for the sake of completeness. (As 
indicated in ref. 3, the lift distribution for a wing of very low aspect 
ratio is independent of the plan form, provided the trailing edge is not 
reentrant. ) The lift distributions on the nineteen wings considered in 
the present paper are presented in figures 2 to 20. Table I serves as 
a table of contents for this group of figures. 
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Aerodynamic Parameters 


The aerodynamic parameters C 


V 


■'EM* 


y*, c 


V 


C 7 , and 
*d' 


L l /2 


for the nineteen plan forms considered are compiled in table II. The 
values of and for a unit effective flap deflection are pre- 
sented in table III, and the values of ^L-jyg ^ or a 

effective aileron deflection are presented in table IV. 


Aerodynamic Influence Coefficients 
for Stipulated Stations 


Aerodynamic influence coefficients for symmetric and antisymmetric 
lift distributions were obtained as shown in appendix A and are pre- 
sented as the matrices an< l [^a] •'• n tables V(a) and V(b), respec- 

tively. Each influence -coefficient matrix in the table applies to a 
given plan form. These influence -coefficient matrices are used to cal- 
culate the spanwise lift distribution for any continuous angle -of -attack 
condition from the following matrix expressions: 


and 


{**■} ■ {“*} 
{ r *a} = C Zd [Qa] {a a } 


( 1 ) 


for the symmetrical and antisymmetrical distributions, respectively, 
where a is the angle of attack at stations y* = O. 9808 , 0 . 9239 * 
0.3315, 0.7071, 0.5556, O. 3827 , 0.1951, and 0 and V* is the desired 
lift at these stations. In this paper the convention is that the angle 
of attack for the station nearest the wing tip (y* = 0 . 9808 ) is the 

first element of the angle-of-attack matrix and the lift at the 

same station is the first element of the lift -distribution matrix {r*} 

The matrices fQgj and JqJ of table V are arranged accordingly. 


DEVELOPMENT OF AERODYNAMIC INFLUENCE COEFFICIENTS 
FOR ARBITRARY STATIONS 


The influence coefficients described in the preceding section are 
satisfactory for many purposes; for instance, the stipulated stations at 
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which the lift is given are convenient for. plotting spanwise lift distri- 
butions because the points are concentrated near the wing tip where the 
curvature of the lift distributions is greatest. In some cases, however, 
other considerations may determine the points on the span at which the 
lift is to be calculated. For instance, when the influence coefficients 
are to be used in an aeroelastic analysis, the location of the stations 
may be dictated by the structural characteristics of the wing} also, if 
lift distributions are to be calculated for the sake of comparison with 
experimental results, this comparison can be facilitated by calculating 
the lift at the same stations at which it is measured and thus avoiding 
the necessity of graphical or numerical interpolation. 

The following sections describe several methods for developing 
aerodynamic influence coefficients for arbitrary stations. 


Method Using Interpolating Matrices 

One way of constructing an influence-coefficient matrix for any 
stations from the matrices presented herein is to calculate interpo- 
lating matrices which give the angles of attack at y* = 0.9808, 0. 9239, 
O.8315, O.707I, O.5556, O.3827, 0.1951* and 0 in terms of the angles of 
attack at the given stations and the values of the lift at the given 
stations in terms of those at the stations y* = O.9808, 0.9239* O.8315* 
O.707I, O.5556, O.3827, 0.1951, and 0. The desired influence-coefficient 
matrix would then be obtained by postmultiplying the one given herein by 
the angle -of -attack interpolating matrix and premultiplying it by the 
lift interpolating matrix. 

In order to illustrate the nature of these calculations, a pair of 
matrices are calculated for stations at every tenth of the semispan by 

linear interpolation. (The stations at which the angle of attack is 

given and those at which the lift is to be found need not be the same, 
but here, as in most cases, they are chosen to be the same as a matter 
of convenience. ) 

With linear interpolation, 

ao.9808 = °- 808 a 1>0 + 0.192OQ.9 

a 0.9239 = 0,2 59^i.o + °’^ la o.g 

a 0.8315 = 0 - 515a O. 9 + °*685Qo.8 


Oq = Oq 
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which can be written in matrix form as 


c 

“0.9808 


0.808 

0.192 

0 

0 

0 

0 

0 

0 

0 

0 

0 


“1.0 

“0.9259 


0.239 

0.761 

0 

0 

0 

0 

0 

0 

0 

0 

0 


“ 0.9 

“0.8515 


0 

0.515 0.685 

0 

0 

0 

0 

0 

0 

0 

0 


“ 0.8 

“0.7071 

< 

> = 

0 

0 

0.071 

0.929 

0 

0 

0 

0 

0 

0 

0 


“0.7 

“0.5556 


0 

0 

0 

0 

0.556' 

0.444 

0 

0 

0 

0 

0 


“0.6 

“0.5827 


0 

0 

0 

0 

0 

0 

0.827 

0.173 

0 

0 

0 

< 

1 “0.5 * 

“0.1951 


0 

0 

0 

0 

0 

0 

0 

0 

0.951 

0.049 

0 


Q 

0 

V- 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.000 


“0.5 


“o.i 


where the rectangular matrix on the right side of the equation is the 
desired angle -of -attack interpolating matrix. 

Similarly, for the lift distribution, 

r*i.o = o 

(This information is known from physical considerations; to calculate 
r*1.0 extrapolation from r* 0 , 9808 and ^*0. 9239 would give a 
spurious value. ) 

r 0.9 = r 0.8315 + 0.0924 V 0-9239 " ^0.8315) 

= o. 74 ir* 0 5239 + °- 2 59r* 0# 83i5 

r 0.8 = r 0.7071 + 0il2 ^V r 0.8315 ‘ r 0.7071J 

= o.747r* 0> 83i5 - 0. 253 r* 0 . 7071 



* 

0 
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or, in matrix form, 


r 

r* 

1 1.0 


0 

0 

0 

0 

0 

0 

0 

0 

* 

r 0.9 


0 

0.7^1 

0.259 

0 

0 

0 

0 

0 

CO 

d 

*L 


0 

0 

0.7^7 

0.253 

0 

0 

0 

0 

1*0.7 


0 

0 

0 

0.953 

0.047 

0 

0 

0 

r*o.6 


0 

0 

0 

0.293 

0.707 

0 

0 

0 

< r *0-5 

S = 

0 

0 

0 

0 

0.679 

0.321 

0 

0 

r *0.4 


0 

0 

0 

0 

0.100 

0.900 

0 

0 

1*0.3 


0 

0 

0 

0 

0 

0.559 

0.441 

0 

r 0.2 


0 

0 

0 

0 

0 

0.026 

0.974 

0 

r* 

r 0.1 


0 

0 

0 

0 

0 

0 

0.513 

0.487 

F* 

1 0 

i 

0 

0 

0 

0 

0 

0 

0 

1.000 






r *0.98o8 

r* 

r 0.9239 

1*0.8315 

r * 

r 0.7071 
r*o.5556 

1*0.3827 

* 

r 0.1951 




l r *° 


I 


where the rectangular matrix on the right side of the equation is the 
desired lift interpolating matrix. 


Although linear interpolation is by far the simplest type, the 
results obtained with it are not so accurate as a higher-order interpo- 
lation procedure. Parabolic interpolation should be satisfactory for 
the angle -of -at tack matrix and for most of the lift-distribution matrix, 
except near the wing tip. The numerical factors required for parabolic 
interpolation can be calculated by means of Lagrange's interpolation 
formula for an mth degree polynomial 


rn-t-1 

*1 n ' (A - A) 
f(y*h) - x f(v\)r 


k=l 


mfl 

n' 

i=l 


(A - A) 
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where f(y*) represents either the lift or the angle -of -attack distri- 
bution, y*^ is the station at which f (y*) is to be determined by 

interpolation, and y*^, y* 2 > • • • are the stations at which f(y*) 

is presumed to be known. The prime mark on the product signs is the 
conventional designation of the fact that the term for i = k is to be 
omitted. For parabolic interpolation (m = 2 ) this formula reduces to 


f(y* h ) = f(y» 1 ) 


to) 


(y^ h - y* 2 )(y*h - ^*3) 


to) 


(y*i - y * 2 ) 
(y* h - y*i) 

( y *i - to 
to - to 

to - to 

to - y *i) 

to - y*,) 

to - to 

'to ■ to 

(y*, - y* 2 ) 


Near the wing tip the lift distribution cannot be approximated 
accurately by an ordinary parabola but can be represented instead by a 

linear superposition of the two functions (l - y *) 1 ^ 2 and (l - y * )^ 2 
as suggested by V. M. Falkner. With this approximation the desired inter- 
polating factors for y* between 0-9239 and 1 are the two elements of the 
row matrix obtained by postmultiplying the row matrix by the square matrix 
in the equation 


r*(j*h) - 1 i 1 - to 1/2 ( x - A ?' 2 


9.652 -1.223 
.126.835 63.708 




r *o. 9808 


0.9239 


For y* h = O. 96 , for instance, the factors obtained in this manner are 
0.9157 and O. 265 I, so that 

4*0. 96 “ 0- 9157^0.9808 + 2 65ir*o.9239 

Lagrange’s general interpolation formula can be used for higher- 
order interpolation (m > 2 ), but the effort entailed in calculating the 
interpolating factors is not generally justified by the increase in 
accuracy obtainable compared to parabolic interpolation. 
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Interpolating factors for a and T* can also "be obtained by 
representing these functions by Fourier series in -fl. For the lift dis- 
tribution the trigonometric interpolation formula 

r *K) = ~4 t 2— f (*i) ^ sin jl8 i sin 

n + 1 i=l j=l 

can be used, where n = 15 for the calculations in this paper, and 
where, as in this paper, the values of are at equal increments 

d, = — . For the angle -of -attack distribution this formula is not 

x n + 1 

applicable, because the angle -of -attack distribution cannot be repre- 
sented accurately by a finite sine series and because the angle-of - 
attack values are presumed to be given at nonequal increments in 
(if they were given at eight equal increments no interpolation would be 
required. ) Although in principle a matrix-inversion method could be 
based on an expansion of the angle-of -attack distribution in a cosine 
series, the matrices to be inverted are generally ill-behaved so that 
the results are likely to be of doubtful accuracy. 


Methods Using Green's Function 

Basis of methods .- Although the aerodynamic -influence -coefficient 
matrices discussed in the preceding section have the property that when 

postmultiplied by the angle -of -attack matrix and multiplied by the 

lift-curve slope they yield the lift-distribution matrix their 

individual elements have no direct physical significance. On the other 
hand, a structural influence coefficient has the significance that it 
represents the deformation at one point caused by a unit concentrated 
force at another point. A corresponding type of aerodynamic influence 
coefficient would represent the lift at one point y* due to a "unit 
concentrated angle of attack" at another point y Q *. This angle of 

attack is actually an angle-of -attack distribution represented by a 
Dirac delta (impulse) function of the distance along the span, that is, 
a function which, in the limiting case as Ay 0 * approaches zero, is 

zero everywhere except in the interval y Q * = y* ^ y Q * + £y 0 * where 
the ordinate is equal to — - — ; the area under this function would 

Ay 0 * 

always be 1, which justifies the use of the term "unit" in connection 
with this distribution. The desired influence coefficients would then 
be the values (at various values of y*) of the lift distributions due 
to these angle -of -attack distributions for various values of y 0 *. 
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The desired lift distributions, which constitute a type of Green's 
function for Weissinger's integral equation, can be obtained as follows. 
Let r*(y*,y 0 *) be the lift distribution for a unit effective deflec- 
tion of a flap which is located between y* = y Q * and y* = 1. Then 

the lift distribution for an angle -of -attack distribution for which 
the angle of attack is zero everywhere except in the interval 

y Q * ^ y* = y 0 * + Ay 0 *> Where it is — *s given by 

r*(y*,y 0 *) - r*(y*,y 0 *+£y 0 *) 

Ay 0 * 


The lift distribution corresponding to the unit concentrated angle-of- 
attack distribution therefore is the limit of this expression as Ay 0 * 

dr*(y*,y 0 *) 

approaches zero, which is by definition. For any given 

^o* 

angle -of -at tack distribution the lift distribution can be determined by 
linear superposition of lift distributions of the Green's function type 
as follows: 


r*(y*) = - 



dr*(y*,y 0 *) 


d y 0 * 


( 2 ) 


The desired Green's function can thus be obtained by calculating 
r*(y*,y 0 *) and taking its partial derivative with respect to y Q *. A 

more convenient approach, however, is to consider symmetric and anti- 
symmetric loadings separately. By a repetition of the preceding argu- 
ment the following results are then obtained: 


r * s (y*) 



(3) 



Q^o*) 


ft*aii(y**y 0 *) 


d y Q * 


(4) 


where, as elsewhere in this paper, r*^ is the lift distribution (as a 
function of y*) for an inboard flap extending over the interval 
_y o * g y* ^ y Q *, and r* a: Q is the lift distribution for outboard 
ailerons extending over the intervals y Q * ^ |y*| ^ 1* The argument 
y Q * in a(y 0 *) can be regarded simply as a variable of integration 
corresponding to y*. 
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The desired Green's function can thus "be obtained from the flap and 
aileron distributions given herein by differentiation with respect to 
y Q *. The results presented in parts (c) and (d) of figures 1 to 20 can, 

for instance, be cross -plotted as functions of y Q * ^y 0 * - — for sym " 

metrical loadings, and y Q * = 1 - for anti symmetrical loadings^ for 

given values of y* and the differentiation then performed graphically. 
The graphical procedure, however, is tedious and relatively inaccurate. 
Similarly, numerical differentiation of these lift distributions would 
tend to be inaccurate. Two numerical procedures which avoid differenti- 
ation of the lift distributions are therefore presented in the following 
sections; one consists in calculating the desired Green's functions 
directly, and the other consists in using derivatives of the angle-of- 
attack distribution. 


Direct calculation of Green's functions . - Inasmuch as the desired 
Green's functions are lift distributions corresponding to angle -of - 
attack distributions defined by delta functions, they can be calculated 
directly provided the singularities in the angle -of -attack and lift dis 
tributions are taken into account. Appendix B of this paper describes 
the method by which the singularities in the flap- and aileron-type 
angle -of -attack and lift distributions were taken into account. This 
method is extended to the case of impulse -type angle -of -attack distri- 
butions in appendix C. 


The resulting lift distributions 

cr* f (y*,y 0 *) 

with the Green's functions 

^o* 


rV 


and - 


and r* a ' are identical 

***n(r*’ y o*) 




They 


contain logarithmic singularities at y 0 * = y* which must be split off 

before the integrals in equations (3) and (4) can be evaluated numeri- 
cally. For a symmetric distribution the procedure is as follows. For 
a given value of y*, 


r *s (y* ) 



a s (y*) f 1 rV(y*,y 0 *)ay 0 * (5) 

J 0 

The first integral does not contain any singularity; the integrand is 
zero at y Q * ~ y*- secon< ^ integral can be evaluated explicitly and 

a s (y* )r* Sl (y*) 
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where r* s ^ (y* ) is the lift distribution for a = 1 over the entire 

span, so that the right side of equation ( 5 ) can be evaluated numerically 
without difficulty by using any set of integrating factors appropriate to 
the stations of interest, such as those of Simpson's rule if the points 
are equally spaced and the number of intervals between y* = 0 and 
y* = 1 is even. 


If these integrating factors are written in the form of a diagonal 
matrix and designated by []f|, equation ( 5 ) can be written as 

r * s (y* ) = [^B'(^yo*^l/l{°B(yo*) - a s (y*)j* + a s (y*) r * Sl (y*) ( 6 ) 


Now, if the row matrix 



is defined as 


|ly*J = [00000100 



with the element 1 at the position corresponding to y* and zeros 
elsewhere, if the matrix j^y^] defined as a square matrix all the 

rows of which are equal to ^Ly*J, and if |_1~] represents the unit 
matrix, then 

r* s (y*) = |_ r *s '(y*'yo*)J 1/1 [bn - [V]] { a s/o*)} + 

r* S;L (y* )[V*j{ a s(y 0 *)] > 

[yj + t*s '(y*^o*)J l/l l/l - [v] { a s(y 0 *)} 


or 


r*s(y*) - 


(7) 


where the matrix |^QgJ defined hy 

y* 


l5d 


y* 


r* Sl (y*) 

% 


lyj 


r *s '(y*>y 0 * ) 

% 


tl[hl - [V]] 


( 8 ) 


i 
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is the row matrix corresponding to y* of the desired influenee- 

r r* s (y*) 

coefficient matrix Qg . The values of can he obtained 

directly from the curves labeled "Constant" of parts (a) or the curves 

for — = 1.0 of parts (c) of figures 1 to 20, and the calculation of 
b 


rY(y*’ y o‘) . 


'h. 


is described, in appendix C. This calculation must be 


repeated for all other values of y* to obtain Q s , so that, finally 


M- 


(y*) 


'let 


r s (y*>y Q *) 




hi 


I. 


(9) 


where 


z 


is a diagonal matrix in which the elements are the sums of 


the elements in the rows of the matrix 


r* s '(y*,yo*) 

% 


1/1 


Similarly, for an antisymmetric distribution the singularity can be 
split off in several ways, one of them being the following. For a given 
value of y* equation (4) can be written as 


r* a (y*) = 


1 (a a (y 0 *) - <^(y*^ r Y(y^y 0 *)dy 0 * - a a (y*) r *a/y*) 


where r* ai (y*) is the lift distribution for a unit effective displace- 
ment of a full-span aileron. In matrix notation this relation may be 
written as 


r *a(y*> = Cj d |Q a J 


y* 



( 10 ) 
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where 


LSI 


y* 


The values of 




r* ai (y*) 


H 


r *aV-yo*) 


L*1 


L 1 ! - fv 


can be obtained from the curves for 


J ail 




( 11 ) 


= 1.0 


of parts (d) of figures 1 to 20. This procedure must again be repeated 
for each value of y* to obtain all the rows of j^Q^j , with the result 


that 


HlI = - 

r\(y*) 


r*aV>y 0 *) 

III - 


l>J 

O 

1 

T 

% J 

L 1 1 

L^a 


( 12 ) 


where the diagonal matrix 
of the rows of the matrix 



consists of the sums of the elements 


r»a (y*,y Q *) 


Lil 


In order to indicate the extent of the effort involved in calcu- 
lating these influence-coefficient matrices, a step-by-step summary of 


these calculations is given for 
this procedure also applies to 


M- 


with the obvious modifications 
In the first six steps Green's 


functions are calculated in accordance with the procedure indicated in 
appendix C; in the remaining steps j^Qgj is calculated in accordance 

with the previous discussion in this section. 


(l) For the values of 0 Q = cos“l y Q * of interest, the values of 
sin n0 o (for n = 1, 3> 5 > • • • 15 ) are obtained from trigonometric 
tables and assembled in a matrix {sin n0 o "l. 
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(2 ) This matrix is premultiplied by the matrix Jjcos nf] . 

(See table VI. ) 

(3) The resulting matrix is premultiplied by the matrix [Fa] con- 

taining the elements F(r)*,y*) defined in appendix A. The evaluation 
of this [ja] matrix is probably the most time-consuming part of the 

calculation because it does not lend itself very readily to high-speed 
automatic computation. For any one of the nineteen plan forms considered 

in this paper, the matrices 2 [pJ and 2 [fJ] are available upon request 
from the National Advisory Committee for Aeronautics. 

(4) The resulting matrix is premultiplied by the matrix [Qs] given 
in the present paper for the plan form of concern, and the matrix obtained 
in this manner is multiplied by the constant l/rc . 

(5) The values of T* s ' ^ are calculated from equation (C3) for the 


D 


given , values of 0 O and for 


= 22. n = 1 2 

16 


and assembled in a matrix 


r* 1 
1 s 


D 


'hx 


8; they are then 


the columns of which 


divided by 
pertain to given values of 0 Q . 

(6) The matrices obtained in steps (4) and (5) are added to obtain 
the matrix 


r* ’ 
1 s 


(7) Values of 


r* Sl (y*) 


% 


are obtained by reading the values of 


cc 


_ at y* defined by y* = cos — , n = 1, 2, ... 8, from figures 2 

16 

' p 

to 20 for a constant value of a and multiplying them by — ; these 


values are assembled in a diagonal matrix 


r* Sl (y*) 


'hx 
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( 8 ) The matrix obtained in step ( 6 ) is postmultiplied by the 
diagonal matrix of integrating factors [_I~|. 


(9) The diagonal matrix 


I 


is assembled from elements calcu- 


lated by adding the elements in a given row of the matrix Obtained in 
step ( 3 ). 

(10) The matrices obtained in steps ( 7 ) and ( 8 ) are added to each 
other and that obtained in step ( 9 ) is subtracted from the sum. The 

resulting matrix is the desired matrix |Qg , as defined in equation ( 9 ). 


The entire calculation thus involves the calculation of 8P values 
sin n0 Q , 

of and of r* s ^ (P being the number of stations y 0 *), 

four matrix multiplications, and three matrix additions, as well as the 
calculation of the 64 elements of the matrix jj’J] , if this matrix is 
not available. 


Calculation using spanwise derivative of the angle -of -at tack distri- 
bution .- Equation ( 3 ) can be integrated by parts to yield 


r* s (y* ) = a s (i)r*f (y*,i) 



fos(yp*) r * 


f(y*,y 0 *) d y 0 


* 


( 13 ) 


where r*f(y*,l) is the lift distribution for a unit effective deflec- 
tion of a full-span flap. Similarly, integrating equation (4) by parts 
yields 

r *a(y*) - £ ^pr* all ( y*,y 0 *)dy 0 * (14 ) 

Here again the argument y 0 * in a(y 0 * ) is merely a variable of inte- 
gration, corresponding to y*. In these equations the lift distri- 
butions and r* a; Q serve as influence functions, so that no Green's 

functions need be calculated. Furthermore, neither these influence 

M y o*) 

functions nor, in most cases of interest, the functions have 

^o* 

singularities in the range of integration, so that the numerical evalu- 
ation of the integrals of equations ( 13 ) and (l4) can be effected very 
readily. 



20 


NACA TN 3014 


With a set of integrating factors [i] for the stations of interest 

pi ^ct a (y 0 *) . . 

I -I • -l J • J / "I \ „ , / r\ \ . / ^ I 17 I rtO-r-4 


and the identity a R (l) = a R (0) + 


'0 


dy Q *, equation (13 ) can he 


written as 


{r* s (y*)} = a s (o){r* f (y*,i)} + [f*f(y*A)] - 




where [ r *f (y*,i)] is a matrix all the columns of which are equal to the 

column matrix , and L 1 ! is the diagonal matrix consisting 

of the integrating factors. Equation (15) can he rewritten in terms of 
a new influence -coefficient matrix j^Q'J defined by 

[«■.] - [r^i - <*> 

l j % % 


{r* s (y*)} - os(o){r* f ( y »,i)]. + 


&<s(y 0 * ) 


Similarly, with a new influence-coefficient matrix a defined 


[»'.] - r<all ^’ Yo - ) ] m 
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equation (l4) can be written as 




( 19 ) 


The matrices ^Q' g J an( * [q'J are * >ase< * on ^he assvunption that 


M y o*) . 


dy * 
J o 


is nonsingular and continuous. However, this assumption is vio- 


lated when a g is discontinuous. Discontinuities in a result from 

control-surface deflection or from deflection of parts of the wing rela- 
tive to the rest of the wing and can he treated in the manner Indicated 
in appendix B or, more simply, by superposition of the lift distributions 
given in figures 1 to 20. In the angle -of -attack distributions for which 
influence coefficients are particularly useful, namely those due to 
structural deformations, discontinuities cannot occur. 


M y o*) 

Discontinuities in can arise if simple beam theory is used 

^o* 

for wings with discontinuous stiffness distributions. Actually this 
theory is inapplicable for such wings, and the spanwise slope of the 
twist is never discontinuous. If, however, simple-beam theory is to be 
used anyway for the sake of convenience and because the errors involved 

are considered to be acceptable, then the matrices £q' s J an< * 

can still be used provided one of the stations is located at the point 
of the discontinuity in the stiffness distribution and provided suitable 
integrating factors are used. 


The objection may be raised against the influence-coefficient mat- 


rices 




that they do not actually express the lift 


distribution in terms of the angle -of -attack distribution but rather 
require its derivative. Inasmuch as the angle -of -attack distribution 
can always be reduced to a continuous one by splitting off the discon- 
tinuous part and treating it as described elsewhere in this report, the 
derivative of the angle -of -attack distribution can be obtained numeri- 
cally by using numerical differentiating factors obtained from any text 
on numerical analysis. These differentiating factors can be assembled 
into a differentiating matrix which when postmult iplied by the matrix 
of the angle -of -attack values yields a matrix of values of the spanwise 

derivative of the angle -of -attack distribution. The matrices [q'sJ 
[?'a] can then be postmultiplied by this differentiating matrix 


and 
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in order to obtain new influence -coefficient matrices which express the 
lift distribution directly in terms of the angle -of -attack distribution. 
However, the main advantage of the method outlined in this section is 
that in aeroelastic calculations, for which aerodynamic influence coef- 
ficients are primarily intended, the angle -of -attack distribution usually 
is obtained by integrating its derivative; the use of the derivative then 
actually saves a calculation. 


In such aeroelastic calculations the lift distributions can be con- 
sidered to consist of a known "rigid wing" part (due to airplane attitude 
or motion, built-in twist, or control deflection), which can be calculated 
initially with due regard to all discontinuities, and an initially unknown 


. part due to structural deformation; the matrices 


[Q's] and 



can 


be used to advantage in calculating the latter part. The calculation of 
a(y*) can then be obviated altogether, because if the structural defor- 
mations are referred to the plane of symmetry the structural part of 
a(y*) is zero for y* = 0, so that the first term on the right sides of 
equations ( 15 ) and ( 17 ) represents a known rigid-wing lift distribution 
and can be included with the others. Thus, in general, 




rigid wing 





structural 


^n* 



The separate treatment of these two parts in an aeroelastic analysis 
presents no difficulties and can be effected in a manner similar to that 
employed in reference 4 for the lift distribution due to aileron deflec- 
tion; the use of [q'^] rather than [qT| requires only the omission of 
one of the integrating matrices in the methods of references 2 and 4. 


DISCUSSION 

General Limitations of the Results of This- Paper 


The Weissinger L-method, its range, and validity have been dis- 
cussed in several previous papers (for example, refs. 3 and 5 ) so that 
in this section only a few comments are made about special applications. 

Number of control points .- In references 1 and 5> four control 
points on the semispan were found to give satisfactory accuracy for the 
lift distributions due to constant angle of attack. In the present 
paper, interest is centered primarily on the lift distributions due to 
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twist and control deflection, and for these cases the additional effort 
entailed in using eight rather than four control points was "believed to 
be warranted by the resulting increase in accuracy. 

Fuselage, nacelle, and tip-tank interference . - The Weissinger 
method and all results presented apply only to wings without fuselages, 
nacelles, or tip tanks. At low angles of attack the lift distribution 
on the wing is not affected to a large extent by the presence of the 
fuselage except when it covers a large part of the wing; the effect is 
largely localized near the wing root and is most pronounced for the 
constant angle -of -attack and flap-deflected cases. For the lift dis- 
tributions due to twist and aileron deflection the presence of the 
fuselage can probably be ignored in most cases. 

Nacelles also tend to affect the lift distribution primarily in 
their own vicinity, but these effects may be significant even for the 
lift distributions due to twist and aileron deflection. Tip tanks tend 
to increase the lift over much of the outer part of the wing to a large 
extent, particularly in the case of lift distributions due to twist and 
aileron deflection. Except for wings with very high or very low aspect 
ratio, these effects can hardly be underestimated and must be taken 
into account in designing the wing. 

High angles of attack . - Potential flow breaks down at high angles 
of attack, and the higher the Mach number the lower the angle of attack 
at which linearized potential -flow theories such as the one employed in 
this paper fail to predict the lift distributions accurately. However, 
the critical design loads often occur at high angles of attack. The 
only suggestion that can be made is that once the rigid -wing lift dis- 
tributions at high angles of attack are estimated on the basis of tests 
or experience, the changes in these distributions due to aeroelastic 
effects can be estimated by means of the results calculated in this 
paper. This procedure cannot be justified theoretically because, 
although the nature of the mutual induction effects between various 
parts of the wing after the flow has separated is still substantially 
the same as before, the lifts caused by these induction effects are not 
those predicted by potential -flow theory. However, the changes due to 
aeroelastic action are small unless the speed is near the flutter or 
divergence speed, so that certain inaccuracies can usually be tolerated 
in estimating them. 

Longitudinal location of the center of pressure .- The Weissinger 
theory yields no information regarding the location of the chordwise 
center of pressure; however, the assumption that at each spanwise sta- 
tion the section center of pressure of the angle -of -attack loading on 
the two-dimensional airfoil section is unchanged in three-dimensional 
flow has been found by lifting-surface calculations (ref. 6, for example) 
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to be largely justified for swept and unswept wings of moderate and- high 
aspect ratio (except near the root and tip). If this assumption is used, 
the longitudinal location of the wing center of pressure may be estimated. 
For low-aspect-ratio wings, the chordwise location of the center of pres- 
sure cannot be determined simply, and lifting-surface methods must be 
used. For flap and aileron deflections, accurate theoretical methods for 
calculating the longitudinal location of the center of pressure are not 
available, but the approximate methods suggested in reference 7 may he 
applied to obtain qualitative information. 

E ffective angle of attack for flap deflection .- In order to deter- 
mine the loading due to flap deflection for wings of high and medium 
aspect ratio (for example, A >4), the effective angle of attack Og 

for the flap (or aileron) deflection may be approximated satisfactorily 
by the values obtained from two-dimensional thin-airfoil theory. Fig- 
ure 21 gives a plot of the effective angle of attack Og against flap- 

chord ratio c f /c. For very low aspect ratios (approaching 0 and cer- 
tainly less than l/2) values of Og close to 1 are indicated by line- 
arized potential-flow theory, even for relatively small values of c f/ c > 

For aspect ratios from about l/2 to 4, lifting-surface methods must be 
used to obtain potential-flow solutions for the lift distributions due 
to partial-chord control deflections. 

Calculation of the roll due to sideslip C^. - The loading for the 

case of full-span ailerons = 1 is the same as the loading on a 

wing with dihedral in yaw or sideslip, because in this case the loading 
on the wing is that due to an angle of attack equal to the product of 
the sideslip angle and dihedral angle on one wing, and the negative of 
that angle of attack on the other wing. The value of Cj^og for the 

case of full-span aileron deflection is therefore equivalent to the 
parameter 

Some other stability derivatives can be deduced similarly from the 
results presented in this paper. 


Relative Merits of the Various types of 
Aerodynamic Influence Coefficients 

In this paper four types of aerodynamic influence coefficients have 
been discussed: 

(l) The influence-coefficient matrices presented in table V which 
are obtained by solving Weissinger 's integral equation by numerical 
methods 
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(2) Influence -coefficient matrices obtained from those of table V 
by multiplying them by interpolating matrices . 

(3) Influence coefficients based on Green's function 

(4) Influence coefficients based on flap- and aileron-type lift 
distributions which express the lift distribution in terms of the span- 
wise derivative of the angle -of -attack distribution rather than the 
distribution itself 

The influence coefficients given in table V apply to the stations 
y* = O. 9808 , 0.9239, O. 8315 , 0.7071, O. 5556 , O. 3827 , 0.1951, and 0. If 
these stations can be used in the calculations in which the influence 
coefficients are to be used, the coefficients given in table V are by 
far the simplest to use because they require no further calculations. 

If these stations cannot conveniently be used, one of the other three 
types of influence coefficients has to be calculated. 

The second type of influence coefficients is based on the first 
and requires a premultiplication of their matrix by a lift interpolating 
matrix and the postmultiplication of that matrix by an angle -of -attack 
interpolating matrix. These interpolating matrices serve to relate the 
lift and angle of attack at the stations of interest to those at the 
stations specified in the preceding paragraph. The interpolating mat- 
rices can be constructed in several ways; parabolic (or possibly cubic) 
interpolation is probably the most satisfactory choice for the angles 
of attack, and for the lift distributions either this type of interpo- 
lation (with a modification at the wing tip) or trigonometric interpo- 
lation should be satisfactory. The interpolating factors do not lend 
themselves readily to automatic computation, but the amount of effort 
involved is relatively small. The two matrix multiplications can then 
be performed readily on automatic computation machines. 

The influence coefficients based on Green's functions are similar 
in concept to the commonly used structural influence coefficients. The 
values of the influence functions r*s ' (y*,y 0 *) and r* a ' (y*,y 0 * ) are 

the only aerodynamic influence coefficients discussed herein which indi- 
vidually have physical significance; the first two types of influence 
coefficients have only a collective physical significance in that they 
yield the values of the lift when matrix -multiplied by the angle -of - 
attack values. However, this individual significance of the coefficients 
based on Green's functions is lost once these coefficients are manipulated 
in the manner indicated in equations ( 9 ) and (12) to obtain aerodynamic 
influence coefficients useful in further computations, and the resulting 
influence coefficients have neither more nor less significance than the 
others. The computation of these coefficients requires a relatively 
large expenditure of effort - four matrix multiplications and three mat- 
rix additions, as well as the computation of many values of f*s ^ or 
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* 1 . 

r a d (by substitution of given values of $ and 0 O in eq. (C3)) 

and a few other minor steps. Despite their conceptual attractiveness, 
these coefficients are therefore practically at a disadvantage compared 
to the other much more simply computed influence coefficients. 

The influence coefficients based on flap- and aileron-type lift 
distributions and on the spanwise derivatives of the angle -of -attack 
distributions are probably the simplest to compute (with the exception 
of those presented in table V)j they require only the reading of the 
values of the lift distributions from the figures of this paper at the 
stations of interest and for the flap and aileron spans corresponding 
to the stations of interest, as. well as the multiplication of the mat- 
rix of these coefficients by a diagonal matrix. (For a symmetric dis- 
tribution, a matrix subtraction is also called for. ) As previously men- 
tioned, the fact that these coefficients express the lift distribution 
in terms of the derivative of the angle -of -attack distribution need not 
be a disadvantage and may actually be an advantage. The decision as to 
whether to use the second or the fourth type of influence coefficients 
(once the decision has been made that the stations implied in the first 
(type are Unsuitable') then becomes' largely a matter of individual prefer- 
ence, guided by decisions in any given case as to the relative conven- 
ience of calculating interpolating factors or reading values from the 
figures in this paper, and of using angle -of -attack distributions or 
their derivatives. 


CONCLUDING REMARKS 


Spanwise lift distributions have been calculated for nineteen 
unswept wings with various aspect ratios and taper ratios and with a 
variety of angle -of -at tack or twist distributions, including aileron 
and flap deflections, by means of Weissinger's method with eight control 
points on the semispan. Also calculated by this method were aerodynamic 
influence coefficients which pertain to a certain definite set of sta- 
tions on the span. Three methods for calculating aerodynamic influence 
coefficients for arbitrary stations have been outlined and their rela- 
tive merits discussed. 

The information presented herein can be used in the analysis of 
untwisted wings or wings with known twist distributions, as well as 
in aeroelastic calculations involving initially unknown twist 
distributions. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , May 5 , 1953 . 
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APPENDIX A 

MATRIX FORMULATION OF THE WEISSINGER METHOD 
Calculation of the Lift Distribution 
and Influence Coefficients 


From two-dimensional thin-airfoil theory it can be shown that, if 
all the vorticity of a plane or parabolically cambered airfoil section 
is concentrated at the quarter-chord line, the downwash angle induced 
at the three-quarter-chord line is equal to the geometric angle of 
attack at the three -quarter -chord line. This circumstance leads to the 
Weiss inger L-method in which the lifting vortex is concentrated at the 
quarter-chord line and the boundary conditions are satisfied at the 
three -quarter -chord line. The Weissinger equation (ref. l) can be 
written as 


a(y*) 


i r 1 d i* + a. r 1 F(n*,y*) 4 ,. 

4it J_ i drj* y* - t)* 8it J_i dr]* 


(Al) 


where 


F(T*,y*) = < 


y* _ n* 


1 + tan Al 


fy* - T)* ' 

\ c */2 . 


fy* _ xj*N 
. c */2 , 


- 1>. (n* > 0) 


r / si 

CVi 

OJ 

+ tan aP* + 

' , (y* - t)*\ 

C\J 

o 

\ c*/2 ) 



- lV + 


(tl* ^ 0) 
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and a is the angle of attack or, more specifically, the streamwise 
slope of the me an -camber surface at the three -quarter -chord line. 


Introduction of the trigonometric variables 0 = cos"-*- y* and 
$ = cos -1 n* into equation (Al) yields 


a(9) = 



dr* (3) . dg 

d3 cos 3 - cos 0 


J_ r r(4 ,e) Sillily 

8rt J o d3 


(A2) 


The solution of equation (A2) is effected in reference 1 by the use 
of trigonometric interpolation and integration formulas. An alternate 
solution based on matrix techniques is presented herein which leads to 
the identical results somewhat more simply and suggests a fairly simple 
setup for routine calculations. 

The function f* is approximated by a finite sine series, as in 
reference 1, so that for any value ■3m 


or 


r *m = XI ®n sin ni3 m 
n 



(A3) 


(A4) 


In the calculations of this paper, the values of n in equation (A4 ) 
have been chosen as n = 1, 3> 5> • • • 15 for the symmetrical loadings 
and n = 2, 4, 6, . . . 14 for the anti symmetrical loadings; values of 
3 m were chosen at 


u m 


mjr 

l6 


with m = 1, 2, 3, . . .8 for the symmetrical loadings and m = 1, 2, 

3> • • . 7 for the antisymmetrical loadings. The use of equal increments 
in 3 is essential in the method of reference 1, but any values of 3 
between 0 and n/2 could have been chosen in the matrix analysis used 
in the present paper. This possibility of using arbitrary stations in 
this matrix version of Weissinger's method is an important advantage in 
that it allows the direct calculation of influence coefficients for 
arbitrary stations. Such a calculation is complicated, however, by the 
fact that some of the matrices which must be inverted tend to be ill- 
behaved when nonequal increments are chosen for 3. 
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d f* ( “6 ) 

The value of — — - required in equation (A2) can be obtained 
from equation (A3 ) as 


na n cos ntf n 


so that the first term of equation (A2) becomes 


, pit 2 I na n cos n3 rr. d3 

— / -a 

1+it J o cos $ " cos e m 


1 V- sin n0 m 

- > na_ 

k n sin 0 m 


1_ r dr *m d0 

1 1+n J q d'fl cos t 3. - i 


1 sin n9 m 


cos 9 m [ 1+ sin 0 




The values of can be expressed in terms of f* m (see eq. (A4 ) ) as 


{an} = [ sin n V] _1 { r *m} 


so that 


where 


1_ r n m d0 

1+n: J q d0 cos 0 - cos 0 m 


- M (d 


L J 4 sin 


The matrices 


sin n0. 


- [n], [sin nGj^j, [sin nO^ -1 , and [bJ 


j^sin 0m _) L L ^ L_ l — — 1 

given in table VII (a) for the symmetrical distributions and table VII (b) 
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for the antisymmetrical distributions. As a result of the orthogonality 

of the sine function, the inverse of the [sin n0 m ^j matrix is the same 

as one-fourth its transpose except for the first and last rows. (See, 
for instance, ref. 8. ) 

The second term of equation (A2) can be integrated numerically by 

dr* 

approximating either F or F by a cosine series. Both approxi- 

dd 

mations will yield identical results, and the latter alternative is 
followed here. Thus let 

. . dr*(d) v — . . 

F(d,9) = y b n (0)cos nd 

dd V 

so that 



where 



is the zeroeth row of the inverse of the matrix 


[cos nif] and consists of the same elements written in the form 

of a diagonal matrix. 
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Similarly, for antisymmetrical distributions 



/ v dr a 

F s(M m )— p M 

cit3 



(A12) 


where 




-1 


The matrices |_dJ , > and 

for symmetrical distributions 
distributions. 



and table VII (b) . for antisymmetrical 


Equation (Al) (or its equivalent, eq. (A2)) can now be expressed 
completely in matrix form. For symmetrical distributions, the equa- 
tion is 


or 




(A13) 


and for antisymmetrical distributions, 



[ G a] | r *a} = |«a} ^ 


It should be noted that the Qb] and (~DJ matrices are invariant 
with plan form and that only the [f] matrices need be computed sepa- 
rately for each plan form; all the matrices are independent of the 
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angle -of -attack conditions. A computing form for the elements of the 
F matrices is given in table VIII; this computing form includes pro- 


vision for calculating the load on swept wings. Sample 2 F„ 

a 


and 2 


matrices are shown in table IX. 



Equations (A13) and (Al4) can be expressed as 


{r*} = [g]-^} (A15) 

so that the elements of the matrices 0-0 constitute, in effect, sets 

of aerodynamic influence coefficients. The influence coefficients pre- 
sented in table V for the plan forms treated in this paper are defined 
as 


and 


so that 


and 




(A16) 


The division by 


and by 



has been performed both to facilitate 


interpolation of the coefficients for unswept wings with plan forms other 
than those considered in this paper and for convenience in aeroelastic 
calculations. Inasmuch as the lift distribution is much less sensitive 
to Mach number than is the over -all magnitude of the lift, an influence- 


coefficient matrix or chosen for the average of the sub- 
sonic Mach number range of interest (that is, for the effective aspect 
ratio A\ll - M^ corresponding to that average Mach number ) will serve 
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for the entire range, 
priate values of 

ref. 9 , for instance, 
effects on and 


provided only that for each Mach number the appro 
and Cj are used in equations (Al6). ^See 

for simple methods of estimating Mach number 

V) 


Calculation of the Aerodynamic Parameters Associated 
With the Lift Distributions 

The values of the lift, induced-drag, bending -moment, and rolling- 
moment coefficients can be obtained conveniently by the use of the inte 
grating matrices derived in this section. 

An integrating matrix for the lift coefficient associated with 
symmetrical loadings can be obtained as follows: 

The lift coefficient can be written as 

c L . f /; r*s 

If, as before, y* = cos 0 and r* s = sin n0 > with n “ 3, 

n 

5, . . . 19 , then 


or 



sin n0 sin 0 d0 



c l * A|^j{r*.} 


where 
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aind If- nor 1 ], is the first row of the matrix j^sin nSj^j ^ given 
in table Vll(a); the matrix j^CjJ is 6 iven in table X. 

Similarly, the integration to obtain the bending -moment coefficient 
for symmetrical loadings can be performed as follows: 

The bending-moment coefficient can be written as 
C BM * £/ 0 O* 


A 

bJ q 


n/2 


a n sin n0 sin 0 cos 0 d0 


n 


or 


A/2 2 2 2 2 

8\3 a l + 5 a 3 " 21 a 5 + 45 *1 ~ 77 & 9 + 


2 2 2 \ 
a T T — , + a 

117 11 165 1 5 221 1 5/ 


C BM = A 


tC BM 


where the matrix 


ll \ — £ 

|_ c bmJ 8 [3 5 


2_ 

21 


_ 2 _ 

45 


2 _ 

77 


117 


2 

165 


221 


[sin n9] 


-1 


is given in table X. 

For antisymmetrical loadings, the half -wing lift coefficient C 
is obtained in the following manner: 

The equation for the half -wing lift coefficient is 


" L l/2 2 , 


r* dy* 


0 
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and, if F 
then 


* a = a n sin n0, with n = 2, k, 6, . 
n 


. Ik as before, 


n Tt/2 

Ct = — / a„ sin n0 sin 0 d0 

L i/a 2J 0 4- n 


• - A (? a 2 " h ^ * » a 6 ' k a 8 * » a l0 - ife a l2 * I95 h*) 


or 


Or / = A 
L l/2 


Lc l 


1/2 




where the matrix 


Lc l 


1/2 



1 

2 3 

4 

5 

6 


_5 

~15 35 

‘63 

99 

“l43 


[sin n0] 


-1 


is given in table X. 

Similarly, the rolling-moment coefficient can be obtained as follows 
The equation for the rolling-moment coefficient is 


or 


Cl 



r* a y* dy* 



a n sin n0 cos 0 sin 0 d0 
n 
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where the matrix 

L ic jJ 3 #L& in h' 1 ], 

is given in table X. 

The induced-drag coefficient C™ can be written as 

i 



where 


or 



dr*s(T]*) dr)* 
dii* y* - tj* 



For symmetrical loadings the 
the values of a^r* s so that 



matrix can be used to integrate 



An integrating matrix to evaluate C D for antisymmetrical distributions 
can be set up similarly. However, in this paper, values of C Di were 
calculated only for additional lift distributions. 
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APPENDIX B 

CALCULATION OF LIFT DISTRIBUTIONS FOR DISCONTINUOUS 
ANGLE -OF -ATTACK CONDITIONS 


The method of solving equation (Al) outlined in appendix A relies 
heavily on numerical integration, as does the method of reference 1. 
Discontinuous angle -of -attack distributions therefore cannot be analyzed 
as accurately as continuous ones can, because discontinuous angle-of- 
attack distributions are known (on the basis of knowledge of the lift dis- 
tributions of wings with very low and very high aspect ratios presented in 

cLP ^ 

ref. 3) to give rise to logarithmic singularities in the function — — 

dr]* 

which occurs in the integrands of both integrals in equation (Al) . Nor 
can a discontinuous angle-of -attack distribution be described adequately 
by a small number of points on the semispan; for instance, with stations 
located as they are for the calculations described in this paper any 
inboard flap terminating at a value of y* greater than O.3827 but less 
than O.5556 would, for a unit effective angle-of -attack distribution, be 
characterized by the angle-of -attack distribution 1, 1, 1, 0, 0, 0, 0, 0 
regardless of the exact location of the end of the flap. 

These difficulties can be overcome by using the results obtained by 
solving equation (Al) for the case of wings of vanishingly small aspect 
ratio (see ref. 3), in which case the second integral vanishes. This 
technique is similar to the one used by Multhopp (ref. 10) in connection 
with the Prandtl lifting-line equation to handle discontinuous angle-of - 
attack distributions. The lift distribution T* is considered to consist 

of a "discontinuous" part r*p which is the solution to equation (Al) if 

the second term is neglected and of a correcting "continuous" part T* c , 
so that 



+ r 


* 

D 


where r* D is defined implicitly by 


1 r n d r*D dd 

kit J 0 dd cos d - cos 9 


(Bl) 


(B2) 
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* 

and F q is the correction that must be added to F p to obtain a 
function T* which satisfies equation (Al) for the given discontinuous 
angle-of -attack distribution. The solution of equation (B2) for 

corresponding to the more common discontinuous angle-of -attack distri- 
butions is given in reference 3; specif ically, for inboard flaps termi- 
nating at y* = y Q * = cos - - 1 - 0 O 


r V*’ 9 °) - 1 


(n - 20 o )sin 3 - (cos 3 - cos 0 o )log € 


sin 


3 + 0Q 


3-0, 


sin 


3 + 0r 


(cos 3 + cos 0 o )log e 


cos 


3-0, 


cos 


(B3) 


and for outboard ailerons with inner ends at y* = y 0 * = cos“l 0 O 


r \il( 9 ' 0 o) ■ 7 


3 + 0, 


(cos 3 - cos 0 o )log e 


sm 


sin 


3-0, 


cos 


2 

3 + 0 O 


(cos 3 + cos 0 o )log e 


3 - 0 f 


cos 


(B4) 


If equations (Bl) and (B2) are substituted into equation (A2 ), the 
result is 


_i f* ar* c (^,e 0 ) 

kn 




0 cos - cos 0 8rc J o 


i r n i \arV*, 0 °) / \ 

- ±- / f( 3,0) — 2L — L d3 = R(e,0 o ) 

Hit 'J n 00 

(B5) 



k-0 
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where the function R(9,0 O ) is defined by 


n r n dr\(d,0 o ) 

R ( 9 ’ 9 °> s rJo F(s ' 9) ^7— 


dO 


or, specifically, for flaps 


Rs(6 ’ 9 ° ) E d 0 F » (9 ’ e) — s — d9 


(B6a) 


and for ailerons 


p n dr* D (d,0 ) 

E ^ e ' e °) s rJo F ° (9 ’ 9) — 17 39 


(B6b) 


Comparison of equation (B5) with equation (A2) indicates that equa- 
tion (B5) may be considered to be the Weissinger equation (eq. (Al) or 
(A2)) for the lift distribution T* c on the given wing (the plan form 

of which determines the function F(d,©)) corresponding to an angle-of- 
attack distribution R(©,0 O ). Inasmuch as R(0,9 O ) is a continuous 

function, as is demonstrated presently, equation (B5) can be solved in 
the manner used for equation (A2). If R(0,0 q ) is being evaluated at 

the stations considered in this paper (9 = 2 3 —r> • • • an( i 

V 16 1 6 16 2/ 

if the eight values of R are listed in a column in the order of 
increasing 0, then premultiplication of this column by the matrix [V] 
given in this report and by the appropriate value of Ct yields the 


desired function 
distribution. 


for the given discontinuous angle -of -attack 


As indicated in equations (B6), the function R(0,9 O ) depends on 

the plan form, which determines F(fl,9), and on the position of the dis- 
continuity in the angle -of -attack distribution, which determines 

dr*(fl,9) 


and, hence, 
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For flaps and ailerons, 


* + e o * + e 0 


D f ^ / sin — cos ‘ 

— — = - (n - 29 0 )cos 4 + sin 4 logg + log — 

dl9 * 4 - e n e 4 - 9, 


D ail 4 

— = - 2 sin 0 O + sin 4 log,, 

d4 * e 


4 + 0, 


4 - 0 , 


- log. 


4 + 0, 


4-0, 


so that the functions — — — may be seen to have logarithmic singulari- 

ot3 

ties. The evaluation of r(0,0 q ) from equations (B6) by numerical 

methods is therefore not a trivial problem. A logarithmic singularity 
Is integrable, however, and F(4,0) is always continuous in 4 and 0 
so that R(0,0 O ) must always be continuous. The integration can thus 

be effected readily by expanding F(4,0) in a finite Fourier series as 
follows: 


F a(*> 0 ) =Z_ P a (6)cos n4 (n = 1, 3, 5, . . . 15). 
n n 


F s (4,0) = > P s (0)cos n4 (n = 0, 2, 4, . . . l6) 


Substitution of these expressions for F and those of equations (B7) 

Sr* 

for — into equations (b 6) yields 


R s (e,8o) * ~ II %(«)[(" - 29 0 )l n ♦ J„ + y 


R a(Mo) “ „ (0 H J n - *n) 

2n n 
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where 


I n = / cos d cos nd dd 


'0 


sin 


/ sin d cos nd logg 


0 


K n s r si 

J 0 


sin 


cos 


d + e c 

~ 

1 - e c 
2 

d + e„ 


dd 


> 


sin d cos nd log e 


0 - 9r 


dd 


cos 


These integrals can all be evaluated explicitly and are 


I 0 = o 


T — « 
X 1 " 2 


I n - 0 


(n = 2, 3, 


Jq = it sin 0 O 


J ± = £ sin 20 o 


J = — 
n 2 


sin(n + l)0 o sin(n - l)0 o 


n + 1 


n - 1 


(n = 2, 3, . . . ) 


and 


Kq = -it sin 0 O 

K 1 - J « in 2 0 o 


K„ - (-1) 


n+1 « 


sin(n + l)9 Q sin(n - 1)0^ 


n+1 


n - 1 


(n = 2, 3, • • • ) 


(BIO) 


>(B11) 

) 

_ > 

\ 

► (B12) 


> (B13) 
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With these values for I n , J n , and equations (B9) can he simplified 


R s( 9 ' e °) = ^X Pa n( 0,0 °) Hn ( 0 °^ (“ = 5, 5, • • • 15) 

n 

R a( 9 ' 0 o) = (n= 0, 2, 4, . . . 16 ) 


(B14) 


where 


Hq = 0 


it - 29 0 + sin 29 c 


sin(n + l)0 o sin(n - l)9 ( 


n + 1 


n - 1 


(n = 2, 3, • • •) 


(B15) 


The function r( 0,9 o ) can also he expressed in matrix form by writing 
equations (B8) in matrix form as 


[_F(d,0)J = |p n (0,0 o )J[cos nfl] 


so that 


p n( e > e o) = [F(d,0)J [[cos n{] 


(Bl6) 


where [cos n-3j is the transpose of the matrix [jcos n-fl] . The following 

expressions are then obtained by combining equation (Bl6) with the matrix 
equivalent of equations (Bl4): 


R s( e > e o) = ^f| p a(^ 0 )J[!- cos “*]_'] 1 { H n( 0 o)j 

R a( 0 > 0 o) = ^ [ F s (*» 0 )J |fcoa nfl ] [\ )j- 


(B17) 
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In essence the procedure used in this appendix for calculating lift 
distributions for discontinuous angle -of -attack distributions can thus 
be seen to consist in avoiding the necessity of integrating numerically 
an initially unknown singular function, as required in equation (Al) or 
(A2), by performing, in effect, the integration of the singular part of 
the function analytically (by solving eq. (B2), as in ref. 3) and 
treating only the continuous part of the function numerically (by 
solving eq. (B5)). As part of this procedure, a singular function has 
to be integrated numerically in order to evaluate R(9,0 o )j however, 

this function is known initially because i't is the product of a known 
singular but integrable function and a known continuous function 

so that, by expanding the continuous 

function in a Fourier series, the numerical integration can be effected 
without difficulty. 


dr* 


D 


and F (3,9), respectively 
d3 / 


The values of 



for flap-span ratios of 0.1, 


0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9 and for aileron-span ratios 
of 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0 are given in 
table XI. 
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APPENDIX C 

CALCULATION OF GREEN'S FUNCTIONS 


As pointed out in the body of this paper (see eqs. ( 3 ) and (4)) 
the desired Green's functions are 

dr* f (y*,y 0 *) dr* ail (y*,y 0 *) 

$y 0 * ' $y 0 * 

^ ~X' 

The functions r ^ and T can he obtained in the manner indicated 

in appendix B. However, in order to calculate the desired derivatives 
of these functions, a numerical differentiation would have to be per- 
formed with respect to y Q *. Such a numerical differentiation is inher 

* * 

ently inaccurate inasmuch as P ^ and r .q have singularities. The 
desired Green's functions are therefore best calculated without using 
the calculated values of r*p and P* a ±± and by using, instead, a 
modification of the method of appendix B. 


Let 


rV(y*,y 0 *) 

r *a’(y*,y 0 *) 


dr *f(y*,y 0 *) 

^o* 

fr*aii(r**y 0 *) 


* t -X- 1 

Then T g and T a are. lift distributions corresponding to 
type angle -of -attack distributions of the following forms: 


(Cl) 


impulse - 



Symmetrical loading 


a 



/ 

/ 

/ 

/ 

/ 

✓ 

/ 

/ 



/ 

/ 

/ 

/ 

/ 

✓ 

/ 

✓ 





1 


a 


Antisymmetrical loading 
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Hence, they must satisfy equation (Al) or its equivalent, equation (A2), 
for these angle -of -at tack distributions. 

Again, as in appendix B, the lift distributions can be considered 
to consist of a discontinuous part, which satisfies equation (Al) for 
the given angle-of-attack distribution if the second integral on the 
right side is disregarded, and of a correction part; that is, 


r* ' = r* ' + r* ' 

8 8 d 8 c 


r * ' _ r* ' + r * ' 

a a d a c 


(C2) 


The functions r* s D and T* a D must thus satisfy equation (B2) for 

the given angle-of-attack distributions. By virtue of the linearity of 
equation (B2) and by virtue of the definitions of the functions in terms 

of r* f and r* ail , respectively, r* g ' and r *a' D can be obtained 

by differentiating with respect to y 0 * the solutions of equation (B2) 

given in equations (B3 ) and (B4) for the flap and aileron angle-of- 
attack conditions. Thus 


* ' k sin -fl + sin 0 O 

f s n = r l°S e 

D rt 6in •S - sin 0 O 


4 "I 

r * D ■ « l0g = 


sin(d + 0 O ) 


(C3) 


sin|d - 0 O | 

Similarly, the functions r* s ' and r* ' must satisfy equation (B5), 

c ^ 

where now r( 0,0 q ) is defined for the symmetrical and anti symmetrical 
loadings, respectively, by 

*V D («,e 0 ) 

[ 0 

and 


^( 9 ' e ») ■ if " r * (M> 




d« 




- if' '.(«*•) 


D (^,0 O ) 


dd 


(c4) 
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The evaluation of these integrals can he effected in the manner 
employed for equations (b6) so that 


R s( e < 6 o) = Sr2I P a n (e>H , n( e o) 
n 

R a(«< e o) - ^]T Ps n (e,H 'n( e o) 
n 


(n = 1, 3, 5, 



(n = 0, 2, 4, . . . 1 6 ) 


(C5) 


where 

H 'n( e o) 3 2 sin n0 o 


(C6) 


for all values of n, and where P„ and P e are the same values as 

“n s n 

those used in appendix B. Thus, in matrix form, for given values of 9 
and 0 Q , 



values of y Q * ranging from 0 to 0-9 for symmetrical distributions and 
0.1 to 0.9 for antisymmetrical distributions. 


The desired Green's functions can thus be calculated in the fol- 
lowing way. For a given value of 9 0 , the values of R s (9,9 0 ) and 

^a(®> 0 o) are calculated for eight equal increments of 9 between 0 

and jt/2 from equation (C7 ). These values are then written as 
columns and premultiplied by the matrix of influence coefficients 
tabulated in this paper for the given plan form in order to obtain 


% 


rV c ( e >8 0 ) 


and 



hQ 


MACA'TN 30Hv 


for the given value of 0 O . To these values are added the values of 


rV D (e,e 0 ) r* a ' D (e,e 0 ) 

and — 

% 


* * * 1 

obtained by dividing the values of P s ^ and T a ^ calculated from 
equation (C3) by . and respectively, for this value of 0 O 

and the given values of 0. This procedure yields 


1 s 


'( 9 > e o) 


and 


* 

r a 


'(9, 0 o) 


'hx 


^The division by C-^ and by is performed to facilitate the 

further calculations required to obtain the desired influence coeffi- 
cients, as explained in the body of this paper.) This calculation is 

repeated for all the values of 9 0 for which the Green's functions are 
desired. 
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TABLE I 

INDEX TO SPANWISE-LIFT -DISTRIBUTION FIGURES 


^A.11 wings are unsweptj 


Plan form 

A 

A 

Figure 

— 

Very low 



1 

3H 

1-5 

0 

2 

312 

1.5 

.25 

3 

313 

1.5 

.50 

4 

314 

1.5 

1.00 

5 

315 

1.5 

1.50 

6 

321 

3.0 

0 

7 

322 

3-0 

.25 

8 

323 

3.0 

.50 

9 

324 

3.0 

1.00 

10 

325 

3.0 

1.50 

L 

11 

331 

6.0 

0 

12 

332 

6.0 

• 25 

13 

333 

6.0 

• 50 

14 

334 

6.0 

1.00 

15 

335 

6.0 

1.50 

16 

341 

12.0 

0 

17 

342 

12.0 

• 25 

18 

343 

12.0 

• 50 

19 

344 

12.0 

1.00 

20 
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TABLE II 

FORCE AND MOMENT COEFFICIENTS FOR ANGLE-OF-ATTACK LOADINGS 


Plan form 

Additional loading 

Damping-in-roll loading 

% 

C BM 

y* 

C Di 

Cl d 

Cl 1/2 

511 

i. 8976 

0.7701 

0.4058 

0.7799 

0.1280 

0.4448 

512 

1.9940 

.8551 

.4188 

.8448 

.1405 

.4781 

515 

2.0006 

.8449 

.4225 

.8495 

.1425 

.4855 

514 

1.9782 

.8455 

.4265 

.8505 

.1454 

.4865 

515 

1.9459 

.8555 

.4295 

. 8o4l 

.1458 

.4871 

521 

2.9956 

1.1765 

.5950 

1.0004 

.2126 

.7544 

522 

5.1747 

1.5177 

.4151 

1.0755 

.2501 

.8548 

525 

5.1755 

1.5410 

.4226 

1.0686 

.2584 

.8772 

524 

5.0970 

1.5564 

.4515 

1.0198 

.2642 

.8912 

525 

5 . 0086 

1.5174 

.4579 

.9679 

.2661 

.8948 

551 

4.1171 

1.5557 

• 5774 

• 9950 

.5085 

1. 1248 

552 

4.5581 

1. 7850 

.4110 

I.OO71 

.5895 

1.5555 

555 

4.5205 

1. 8565 

.4250 

.9918 

.4i4l 

1.5964 

554 

4.1816 

1.8479 

.4419 

.9409 

.4554 

1.4417 

555 

4.0252 

1. 8260 

.4556 

.8921 

.4405 

1.4540 

541 

5.0125 

1. 8145 

.5620 

.7850 

.5918 

1.4714 

542 

5.1989 

2.1146 

.4067 

• 7510 

.5190 

1. 7797 

545 

5.1647 

2.2164 

.4291 

.7126 

• 5700 

l. 9008 

544 

5.0026 

2.2819 

.4561 

.6982 

.6156 

2.0024 
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TABLE V 

AERODYNAMIC- INFLUENCE-COEFFICIENT MATRICES 


(a) Symmetric loadings 





Plan form 311 




0.0530 

0.0360 

0.0254 

0.0218 

0.C202 

0.0182 

0.0184 

0.0086 

0.0224 

0.1350 

0.0914 

0.0644 

0.0587 

0.0517 

0.0508 

0.0243 

0.0121 

0.0677 

0.2298 

0.1535 

0.1144 

0.1041 

0.0956 

0.0475 

0.0093 

0.04l6 

0.1268 

0.3275 

0.2221 

0.1722 

0.l6i4 

O.O768 

0.0075 

0.0341 

0.0857 

0.1951 

0.4276 

0.2979 

0.2432 

0.1194 

0.0068 

0.0289 

O.0739 

0.1426 

0.2746 

0.5320 

0.3910 

0.1720 

0.0064 

O.0273 

0.0664 

0.1297 

0.2170 

0.3737 

0.6582 

O.2788 

0.0065 

0.0262 

0.0658 

0.1233 

0.2114 

O.3258 

0.5497 

0.4974 




Plan form 312 




0.0590 

O.O536 

0.0430 

0.0401 

O.O356 

0.0338 

0.0318 

0.0159 

0.0274 

0.1457 

0.1097 

O.O855 

O.O782 

0.0704 

0.0680 

0.0330 

0.0154 

0.0763 

0.2365 

O.167I 

0.1296 

0.1189 

0.1095 

0.0544 

0.0116 

0.0479 

0.1332 

O.3269 

0.2269 

O.I789 

0.1676 

0.0801 

0.0091 

0.0385 

0.0904 

0.1960 

0.4170 

0.2929 

0.2399 

O.U78 

0.0080 

0.0323 

0.0769 

0.1429 

O.2676 

0.5U3 

0.3759 

0.1655 

0.0073 

0.0301 

0 . 0686 

0.1289 

0.2106 

0.3581 

0.6274 

O.2656 

_o . 0072 

0.0291 

0.0676 

0.1223 

0.2047 

0.3118 

0.5233 

0.4733_ 




Plan form 313 




0.0591 

0.0544 

0.0446 

0.0424 

0.0382 

O.0363 

0.0343 

0.01 7 r 

0.0277 

0.1468 

0.1123 

O.O893 

0.0826 

0.0747 

0.0721 

0.0350 

0.0157 

0.0776 

0.2388 

0.1709 

0.1342 

0.1236 

o.n4o 

0.0566 

0.0119 

0.0489 

O.1350 

0.3292 

0.2300 

0.1821 

0.1706 

0.0815 

0.0093 

0.0391 

0.0914 

0.1971 

0.4174 

0.2933 

0.2400 

O.H78 

0.0081 

0.0325 

0.0770 

0.1426 

0.2662 

O.5083 

0.3725 

0.1637 

0.0073 

0.0300 

0.0680 

0.1276 

0.2079 

0.3536 

0.6209 

0.2622 

,0.0073 

0.0288 

0 . 0668 

0.1206 

0.2014 

0.3066 

0.5161 

0.4688 




Plan form Jlk 




0.0599 

0.0554 

0.0460 

0.0442 

0.0401 

0.0382 

0.0361 

O.OI79 

0.0283 . 

0.1492 

0. 1151 

0.0925 

0.0860 

0.0780 

0.0752 

0.0365 

0.0162 

0.0793 

0.2432 

0.1752 

0.1384 

0.1274 

0.1174 

0.0582 

0.0122 

0.0501 

0.1377 

0.3344 

0.2341 

O.I853 

0.1751 

0.0826 

0.0094 

0.0396 

0.0925 

0.1991 

0.4213 

0.2948 

0.2400 

O.H76 

0.0081 

O.0323 

0.0767 

0.1419 

0.2654 

0.5083 

0.3695 

O.1616 

0.0072 

0.0294 

0.0666 

0.1248 

0.2035 

0.3480 

0.6160 

O.2588 

0.0070 

0.0280 

0.0648 

0.1169 

0.1955 

0.2986 

0.5077 

0.4664 
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TABLE V.- Continued 

AERODYNAMIC- INFLUENCE-COEFFICIENT MATRICES 


(a) Continued 


Plan form 315 , 

0.0565 0.0470 0.0452 o.o4io 0.0390 0.0367 0.0182 

0.1519 0.1174 0.0945 O.O878 0.0795 0.0764 0.0370 

0.0808 0.2475 0.1784 0 . 1408 0.1293 O.II87 O.O588 

0.0509 0.1399 0.3398 0.2374 O.187I 0.1740 0.0828 

0.0400 0.0933 0.2011 0.4261 0.2964 0.2397 0.1171 

0.0322 0.0764 0.1415 0.2653 0.5104 O.3677 0.1600 

0.0288 0.0651 0.1224 0.1999 0.3439 0.6138 0.2564 

0.0270 0.0627 0.1132 0.1900 0.2912 0.5005 0.4657 

Plan form 321 

0.0272 0.0128 O.OO87 0.0062 0.0064 0.0047 0.0057 

0.0094 0.0716 0.0402 0.0241 0.0209 0.0174 0.0169 

0.0042 0.0309 0.1245 O.O718 0.0467 0.0405 O.0355 

0.0031 0.0164 O.O607 0.1793 O.IO77 0.0746 0.0675 

0.0023 0.0128 0.0365 0.0961 O.2353 0.1483 0. 111 7 

0.0020 0.0102 0.0302 0.0641 O.I385 0.2945 0.2011 

0.0018 0.0096 0.0262 O.O569 0.1029 0.1944 O.3696 

_0.0019 0.0092 0.0261 0.0534 0.1000 0.1647 0.3014 

Plan form 322 

0.0363 0.0309 0.0219 0.0186 0.0150 

0.0159 0.0864 0.0591 0.o4o4 0.0340 

0.0080 0.04i4 0.1353 0.0855 0.0584 

0.0056 0.0232 0.0688 0.1820 0.1127 

0.0040 0.0173 0.0417 0.0983 0.2281 

0.0034 0.0136 0.0335 0.0648 0.1328 

0.0030 0.0124 0.0285 0.0564 O.O97O 

_0.0029 0.0118 0.0281 0.0525 0.0936 

Plan form 323 

O.O367 0.0323 0.0242 0.0210 0.0171 

O.OI65 O.O887 0.0631 0.0448 0.0379 

0.0086 0.0437 0.1395 0.0907 0.0630 

0.0060 0.0248 0.0720 0.1861 0.1166 

0.0042 0.0182 0.0434 0.1003 0.2296 

0.0035 0.0140 0.0340 0.0650 O.1317 

0.0030 0.0124 0.0284 0.0554 0.0942 

JD.0029 0.0117 0.0277 0.0511 0.0901 



0.0152 0.0136 O.OO67 
0.0314 0.0291 0.0138 
0.0537 0.0465 0.0228 
0.0820 O.O727 0.0337 
0.1444 0.1072 O.0517 
0.2746 0.1831 0.0764 
0.1747 0.3352 O.1357 
0.1448 0.2682 0.2653 


0.0135 0.0122 0.0061 
0.0283 0.0266 0.0126 
0.0501 0.0437 0.0216 
0.0793 0.0709 0.0330 

0.1441 0.1080 O.O523 
0.2774 0.1874 0.0788 
0.1801 0.3427 0.1400 
0.1512 0.2772 0.2703 


0.0022 

0.0080 

O.OI77 

0.0313 

0.0542 

0.0850 

0.1520 

0.2910 


0.0609 

0.0288 

0.0165 

0.0124 

0.0095 

0.0080 

0.0070 

O.OO67 
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TABLE V.- Continued 

AERODYNAMIC-IMFLUENCE-COEFFICIENT MATRICES 


(a) Continued 


0.0379 0.0339 0 
0.0173 0.0922 0 
0.0092 0.0462 0 
0.0064 0.0264 0 

0.0045 0.0191 o 
0.0035 0.0141 0 
0.0029 0.0121 0 
0.0028 0.0111 0 


0.0390 0.0351 0 
0.0179 0.0954 0 

0.0096 0.0480 0 

0.0067 0.0276 0 
0.0046 0.0197 0 
0.0035 0.0142 0 

0.0028 0.0117 0 
0.0026 0.0105 0 


Plan form 324 

.0262 0.0232 0.0190 
.0671 0.0489 0.0415 
.1458 0.0965 0.0678 
.0758 0.1935 0.1218 
.0453 0.1036 O.2351 
.0343 0.0651 0.1314 
.0275 0.0534 0.0900 
.0263 0.0481 0.0843 

Plan form 325 

.0274 0.0244 0.0200 
.0698 0.0513 0.0435 
.1510 0.1005 0.0706 
.0787 0.2000 0.1257 
.0466 0.1063 0.2409 
.0345 0.0652 0.1318 
.0267 0.0517 0.0868 
.0249 0.0453 0.0795 

Plan form 331 


0.0141 

O.OO33 

0.0011 

0.0008 

0.0005 

0.0005 

0.0004 

0.0004 


0.0028 

0.0397 

0.0129 

0.0051 

O.OO38 

0.0027 

0.0026 

0.0024 


0.0028 

0.0155 

0.0709 

0.0273 

0.0129 

0.0101 

0.0081 

0.0082 


0.0010 

O.OO76 

0.0314 

0.1039 

0.0454 

0.0249 

0.0212 

0.0191 


0.0019 

0.006l 

0.0163 

0.0500 

0.1374 

O.O678 

0.0436 

0.0422 


Plan form 352 


0.0255 

0.0098 

0.0040 

0.0024 

0.0015 

0.0012 

0.0010 

0.0010 


0.0189 

0.0565 

0.0228 

0.0102 

0.0068 

0.0047 

0.0041 

0.0038 


0.0108 

0.0322 

O.O837 

0.0351 

0.0170 

0.0124 

0.0096 

0.0095 


0.0055 

0.0132 

0.0235 

0.0548 

0.1340 

0.0642 

0.0395 

0.0377 


0.0080 

0.0175 

0.0432 

O. 1 O 90 

0.0482 

0.0258 

0.0210 

0.0186 


0.0166 0.0147 
0.0341 0.0310 

0.0571 0.0486 
0.0850 0.0740 
0.1460 0.1061 
O.2750 0.1785 
0.1682 0.3283 
0.1352 0.2562 


0.0173 0.0150 
0.0354 0.0316 
0.0588 0.0492 
0.0868 0.6743 
0.1477 0.1050 
0.2772 0.1756 
0.1636 0.3243 
O.1273 0.2461 


0.0006 0.0017 0.0003 
0.0049 0.0045 0.0022 

0.0134 0.0110 0.0055 
0.0285 0.0246 0.0109 
0.0718 0.0462 0.0222 
0.1727 0.1016 0.0394 
0.0993 0.2192 0.0846 
0.0784 0.1693 0.1850_ 


0.0047 0.0040 0.0020 

0.0097 0.0089 0.0041 
0.0184 0.0147 0.0073 
0.0312 0.0261 0.0115 
0.0693 0.0434 0.0206 
0.1607 . 0.0913 0.0346 
0.0884 O.1985 0.0748 
0.0677 0.1493 0.l688_ 



0.0072 

0.0146 

0.0237 

0.0340 

0.0507 

0.0732 

0.1306 

0.2612 


0.0073 

0.0147 

0.0238 

0.0337 

0.0497 

O.0707 

0.1267 

0.2586 
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TABLE V.- Continued 

AERODYNAMIC-INFLUENCE-COEFFICIENT MATRICES 






(a) Continued 








Plan form 355 





0.0263 

0.0211 

0.0135 

0.0100 

O.OO69 

0.0057 

0.0046 

0.0023 


0.0108 

0.0602 

0.0572 

0.0215 

0.0159 

0.0115 

0.0101 

0.0046 


0.0048 

0.0258 

0.0895 

0.0487 

0.0271 

0.0206 

0.0160 

O.OO78 


0.0029 

0.0120 

0.0588 

0.1l4o 

0.0585 

0.0330 

0.0269 

0.0117 


0.0017 

0.0077 

0.0188 

O.O505 

0.1359 

0.0694 

0.0423 

0.0199 


0.0013 

0.0051 

0.0131 

0.0263 

0.0635 

0.1582 

O.O87O 

0.0322 


0.0010 

0.0043 

0.0098 

0.0206 

0.0374 

0.0833 

0.1906 

0.0702 


0.0010 

0.0039 

0.0096 

0.0179 

0.0350 

0.0618 

0.1394 

0.1631 





Plan form 354 





0.0275 

0.0231 

0.0156 

0.0120 

O.OO83 

O.OO65 

0.0051 

0.0025 


0.0118 

0.0644 

0.0420 

0.0256 

0.0188 

0.0131 

0.0110 

0.0048 


0.0055 

0.0289 

0.0968 

0.0550 

0.0312 

0.0229 

0.0169 

O.OO81 


0.0035 

0.0138 

0.0432 

0.1220 

0.0635 

0.0351 

0.0272 

0.0115 


0.0020 

0.0086 

0.0209 

0.0540 

0.l4l4 

O.O705 

0.0407 

O.OI87 


0.0014 

0.0054 

0.0138 

0.0268 

0.0635 

0.1582 

0.0820 

0.0290 


0.0010 

0.0043 

0.0096 

0.0196 

0.0345 

0.0775 

0.1824 

0.0640 


0.0010 

0.0037 

0.0090 

O.OI65 

0.0311 

0.0555 

0.1254 

O.1565 





Plan form 355 





0.0287 

0.0244 

0.0169 

0.0152 

0.0091 

O.OO69 

0.0052 

0.0024 


0.0124 

0.0677 

0.0451 

0.0281 

0.0205 

0.0159 

0.0112 

0.0048 


0.0059 

0.0310 

0.1023 

0.0592 

0.0358 

0.0241 

0.0170. 

0.0080 


0.0036 

0.0150 

0.0463 

0.1286 

O.O672 

0.0565 

0.0271 

0.0111 


0.0021 

0.0092 

0.0222 

0.0567 

0.1468 

0.0717 

0.0394 

0.0177 


0.0014 

0.0055 

0.0140 

0.0271 

0.0657 

0.1597 

O.O786 

0.0265 


0.0010 

0.0041 

0.0092 

0.0188 

0.0524 

0.0752 

0.1777 

0.0594 


0.0009 

0.0035 

0.0083 

0.0149 

0.0281 

0.0475 

0.1146 

0.1520_ 





Plan form 34l 




0.0075 - 

■0.0008 

0.0013 

-0.0003 

0.0008 

-0.0002 

0.0007 

-0.0001 

0.0008 

0.0230 

0.0040 

0.0023 

0.0012 

0.0014 

0.0009 

0.0006 

0.0002 

0.0044 

0.0422 

0.0114 

0.0047 

O.OO36 

0.0029 

0.0014 

0.0002 

0.0012 

0.0109 

O.O631 

0.0207 

0.0089 

' 0.0074 

0.0031 

0, 

.0001 

0.0010- 

0.0035 

0.0196 

0.0845 

0.0319 

0.0159 

0.0077 

0 , 

.0001 

0.0005 

0.0028 

O.OO77 

0.0308 

0.1062 

0.0476 

0.0157 

0.0001 

0.0006 

0.0019 

0.0064 

0.0149 

0.0471 

0.1345 

0.0452 

0.0001 

0.0005 

0.0021 

0.0052 

0.0148 

0.0316 

0.0918 

0.1247 
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TABLE V.- Continued 

AERODYNAMIC-INFLUENCE-COEFFIC IE NT MATRICES 
(a) Concluded 


0.0195 0.0113 0.0046 
0.0059 0.0387 0.0161 
0.0017 0.0115 0.0541 
0.0009 0.0036 O.OI63 
0.0005 0.0022 0.0056 
0.0004 0.0013 0.0038 
0.0003 0.0012 0.0025 
J).0005 0.0010 0.0028 


0.0208 0.0141 0.0066 

0.0072 0.0436 0.0210 
0.0024 0.0146 0.0602 
0.0012 0.0049 0.0195 
0.0006 0.0028 0.0067 
0.0004 0.0015 0.0042 

0.0003 0.0013 0.0026 
_0.0003 0.0011 0.0027 


0.0222 0.0165 0.0088 
0.0084 0.0484 0.0259 
0.0031 0.0178 O.O675 
0.0016 0.0064 0.0233 
0.0007 0.0034 0.0081 
0.0005 0.0017 0.0046 
0.0003 0.0013 0.0027 
0.0003 0.0010 0.0025 


Plan form 342 

0.0029 0.0018 0.0014 

0.0062 0.0043 0.0028 

0.0198 O.OO78 0.0056 
0.0688 0.0246 0.0101 
0.0219 O.O837 0.0310 

O.OO83 0.0290 0.0991 
0.0064 0.0131 0.0403 
0.0051 0.0126 0.0253 

Plan form 343 

0.0042 0.0023 0.0018 
0.0088 0.0057 0.0035 
0.0244 0.0097 O.OO67 
0.0736 0.0273 0.0109 
0.0237 0.0855 0.0311 
O.OO87 0.0286 0.0972 

O.OO63 0.0120 O.O367 
0.0049 0.0113 0.0216 

Plan form 344 

0.0057 0.0031 0.0022 
0.0118 0.0074 0.0041 

0.0298 0.0121 O.OO78 

0.0805 0.0310 0.0119 
0.0262 O.O896 0.0316 

0.0091 0.0284 0.0966 
0.0060 0.0107 0.0323 
0.0043 0.0095 0.0169 


0.0012 0.0006 
0.0025 0.0011 
0.0040 0.0020 
0.0080 0.0032 
0.0145 O.OO69 
0.0413 0.0129 
0.1202 0.0382 
0.0771 0 . 1124 _ 


0.0014 0.0007 
0.0030 0.0013 
0.0044 0.0022 
O.OO83 0.0032 
0.0137 0.0064 
0.0381 0.0113 
0.1135 0.0340 
0.0682 0.1063 


0.0015 0.0007 
0.0033 0.0013 
0.0047 0.0023 

0.0084 0.0031 
0.0126 0.0057 
0.0344 0.0091 
0.1060 0.0284 
0.0557 O.O983 
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TABLE V.- Continued 

AEROEOTAMIC- INFLUENCE -COEFFICIENT MATRICES 


(b) Antisymmetric loadings 



Plan form 311 


0.7812 

0.5201 

0.3395 

0.2619 

0.1930 

0.1255 

0.0642 

0.3225 

1.9623 

1.2583 

O.7818 

0.5867 

O.363O 

0.1869 

0.1615 

0.9296 

3.2298 

1.9474 

1.1733 

O.7708 

O.3634 

0.1102 

0.5007 

1.5966 

4.3377 

2.4478 

1.3241 

0.6582 

0.0703 

0.3351 

0.8640 

2.1287 

5.1126 

2.5737 

1.0511 

0.0447 

0.1977 

0.5300 

1.0690 

2.3425 

5.3269 

2.0671 

0.0210 

0.0969 

0.2420 

0.5089 

0.9126 

1.9539 

4.5926 




Plan form 312 



0.8318 

0.7328 

0.5470 

0.4566 

0.3302 

0.2248 

0.1097 

0.3748 

2.0118 

1.4266 

0.9806 

0.7429 

0.4736 

0.2406 

0.1956 

0.9904 

3.1508 

1.9973 

1.2558 

O.8336 

O.3969 

0.1314 

0.5464 

I.5856 

4.1026 

2.3617 

1.3018 

0.6484 

0.0834 

0.3609 

0.8660 

2.0280 

4.7392 

2.4026 

O.9868 

0.0523 

0.2127 

O.5278 

1.0227 

2.1767 

4.8872 

1.8992 

0.0246 

0.1035 

0.2411 

0.4856 

0.8498 

1.7940 

4.1950 




Plan form 

313 



0.8235 

0.7338 

0.5596 

0.4754 

0.3494 

0.2394 

0.1175 

0.3744 

2.0015 

1.4393 

I.OO87 

0.7730 

0.4976 

0.2534 

0.1974 

0.9935 

3.1400 

2.0136 

1.2827 

0.8564 

0.4096 

0.1328 

0.5509 

1.5877 

4.0796 

2.3641 

1.3117 

0.6545 

0.0841 

0.3629 

O.8678 

2.0191 

4.6940 

2.3853 

0.9817 

0.0525 

0.2130 

0.5268 

1.0168 

2.1546 

4.8256 

1.8756 

0.0245 

0.1032 

0.2399 

0.4814 

0.8400 

1.7702 

4.1376 




Plan form Jlk 



0.8182 

0.7331 

O.565O 

0.4849 

0.3604 

0.2485 

O.1225 

0.3737 

1.9941 

1.4446 

1.0239 

0.7912 

0.5132 

0.2620 

0.1983 

0.9950 

3.1336 

2.0240 

1.3009 

O.8727 

0.4191 

0.1337 

0.5541 

1.5902 

4.0687 

2.3691 

1.3210 

0.6602 

0.0845 

0.3641 

0.8692 

2.0148 

4.6690 

2.3764 

0.9792 

0.0524 

0.2125 

0.5247 

1.0110 

2.1387 

4.7865 

i.8596 

0.0243 

0.1022 

0.2373 

O.476O 

O.83OI 

1.7517 

4.1009 




Plan form 315 



0.8167 

O.7328 

0 . 5666 

0.4880 

0.3640 

0.2516 

0.1243 

0.3735 

1.9920 

1.4465 

1.0289 

0.7974 

0.5186 

'0.2651 

0.1987 

0.9959 

3.1324 

2.0281 

1.3076 

O.8787 

0.4225 

0.1341 

0.5555 

1.5918 

4.0667 

2.3720 

1.3248 

0.6624 

0.0847 

0.3646 

0.8698 

2.0136 

4.6616 

2.3734 

0.9781 

0.0522 

0.2119 

0.5230 

1.0073 

2.1311 

4.7716 

1.8524 

0.0241 

0.1012 

0.2351 

0.4716 

0.8231 

1.7414 

4.0865 
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TABLE V.- Continued * 

AERODYNAMIC - INFLUENCE - COEFFICIENT MATRICES 


(b) Continued 
Plan form 321 


0.3819 

0.1782 

0.1140 

0.0750 

0.0612 

0.0344 

0.0201 

0.1304 

0.9993 

0.5425 

0.2946 

0.2156 

0.1279 

0.0662 

0.0559 

0.4174 

1.7058 

0.9162 

0.4963 

0.3186 

0.1453 

0.0376 

0.2005 

0.7720 

2.3620 

1.2307 

0.6165 

0.3019 

0.0229 

0.1324 

0.3869 

1.0925 

2.8688 

1.3746 

0.5357 

0.0147 

0.0758 

0.2355 

O.5227 

1.2702 

3.0808 

1.1692 

0.0068 

0.0377 

O.IO58 

0.2488 

0.4821 

1.1144 

2.7272 




Plan form 322 



0.4589 

0.3835 

0.2601 

0.2025 

0.1368 

0.0907 

0.0430 

0.1970 

1.0797 

0.7115 

0.4433 

O.3183 

0.1930 

0.0969 

0.0949 

0.4986 

1.6552 

0.9719 

O.5598 

0.3579 

0.1645 

0.0604 

0.2542 

0.7804 

2.1421 

1.1532 

0.5925 

0.2888 

0.0364 

0.l6l4 

0.3982 

1.0016 

2.4972 

1.2055 

0.4721 

0. 0224 

0.0919 

0.2370 

0.4797 

1.1032 

2.6290 

0.9959 

0.0103 

0.0445 

0.1061 

0.2257 

0.4171 

0.9472 

2.3152 




Plan form 323 



0.4488 

0.3874 

O.2771 

0.2216 

O.1525 

0.1008 

0.0479 

0.1978 

1.0708 

0.7326 

0.4758 

0.3453 

0.2108 

0.1053 

O.0983 

O.5072 

1.6481 

0.9954 

0.5873 

0.3760 

0.1733 

0.0628 

0.2628 

0.7886 

2.1165 

1.1555 

0.5994 

0.2917 

O.O376 

O.I654 

0.4033 

0.9926 

2.4374 

1.1795 

0.4622 

0.0228 

0.0930 

0.2370 

0.4725 

1.0718 

2.5423 

0.9604 

0.0104 

0.0445 

0.1050 

0.2200 

0.4019 

0.9107 

2.2341 




Plan form 32k 



0.44l4 

0.3873 

0.2857 

0.2338 

0.1639 

0.1084 

0.0515 

0.1974 

1.0620 

0.7431 

0.4971 

0.3653 

0.2246 

0.1117 

0.1003 

O.5118 

1.6430 

1.0130 

0.6098 

0.3913 

0.1804 

0.0645 

0.2690 

0.7958 

2.1044 

1.1626 

0.6077 

0.2948 

0.0384 

0.1681 

0.4074 

O.9887 

2.4008 

1.1629 

0.4547 

0.0228 

0.0930 

0.2353 

0.4651 

1.0466 

2.4786 

0.9308 

0.0102 

0.0435 

0.1022 

0.2125 

0.3854 

O.8768 

2.1721 




Plan form 325 



0.4388 

O.3872 

0.2886 

0.2384, 

0.1684 

0.1114 

0.0529 

0.1973 

1.0590 

0.7468 

0.5053 

0.3735 

0.2302 

0.1141 

0.1011 

O.5138 

1.6421 

1.0205 

O.6194 

0.3977 

0.1831 

0.0652 

O.2717 

0.7996 

2.1025 

1.1672 

0.6115 

0.2957 

O.O387 

0.1693 

0.4093 

0.9881 

2.3888 

1.1562 

0.4504 

0.0227 

0.0925 

0.2337 

0.4605 

1.0341 

2.4510 

0.9153 

0.0100 

0.0426 

0.0997 

0.2070 

0.3745 

O.8566 

2.1423 
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0.1887 

0.0437 

0.0137 

0.0095 

0.0051 

O.OO35 

0.0015 


.2831 

.1081 

0.0431 

0.0246 

0.0132 

0.0079 

0.0035 


0.2753 

0.1115 

0.0479 

0.0273 

0.0144 

O.OO83 

0.0036 


0.2647 

0.1122 

0.0509 

0.0294 

0.0154 

0.0084 

O.OO35 


0.2615 

0.1124 

0.0521 

0.0303 

0.0158 

0.0084 

0.0034 


TABLE V.- Continued 

AERODYNAMIC- INFLUENCE- COEFFICIENT MATRICES 




(t) Continued 





Plan form 331 



0.0376 

0.0344 

0.0129 

O.OI77 

.0.0057 

O.OO58 

0.5288 

0.2026 

0.0926 

0.0644 

0.0374 

0.0191 

0.1688 

0.9376 

0.3981 

0.1798 

0.1129 

0.0489 

0.0628 

0.3468 

1.3490 

0.5916 

0.2556 

0.1235 

0.0414 

0.1437 

0.5372 

1.7069 

0.7256 

0.2545 

0.0216 

0.0864 

0.2245 

0.6831 

1.9242 

O.6832 

0.0112 

0.0367 

0.1065 

O.2383 

0.6611 

1.8023 



Plan form 332 



0.2092 

0.1158 

0.0809 

0.0486 

0.0316 

0.0142 

0.6255 

O.3496 

0.1787 

0.1188 

0.0661 

0.0331 

0.2476 

0.9187 

0.4544 

0.2184 

0.1327 

0.0571 

0.1046 

0.3690 

1.1751 

0.5408 

0.2382 

0.1129 

0.0615 

0.1585 

0.4750 

1.3884 

0.5908 

0.2060 

0.0321 

0.0899 

0.1971 

0.5470 

1.5138 

O.5296 

0.0155 

0.0378 

0.0907 

0.1867 

0.5093 

1.4097 



Plan form 333 



0.2180 

0.1338 

0.0953 

O.O578 

0.0364 

0.0163 

0.6234 

0.3776 

0.2063 

0.1361 

0.0752 

0.0368 

0.2624 

0.9183 

0.4799 

0.2381 

0.1423 

0.0608 

0.1152 

0.3821 

1.1499 

0.5416 

0.2399 

0.1122 

0.0660 

0.1655 

0.4679 

1.3222 

O.5615 

0.1939 

0.0337 

0.0908 

0.1915 

0.5135 

1.4125 

0.4873 

0.0158 

0.0375 

0.0860 

0.1711 

0.4654 

1.3087 



Plan form 33^ 



0.2201 

0.1451 

0.1065 

0.0657 

0.0403 

O.OI78 

0.6164 

0.3942 

0.2282 

0.1509 

0.0827 

0.0394 

0.2715 

0.9175 

0.5017 

0.2568 

0.1508 

0.0635 

0.1235 

0.3941 

1.1380 

0.5480 

0.2431 

0.1113 

0.0694 

0.1716 

0.4659 

1.2783 

0.5398 

O.1827 

0.0342 

0.0906 

0.1860 

0.4856 

1.3314 

0.4481 

0.0153 

0.0359 

0.0802 

0.1549 

0.4220 

1.2210 


Plan form 335 



0.2207 

0.1494 

0.1114 

0.0692 

0.04l8 

0.0181 

0.6139 

0.4006 

0.2378 

0.1577 

O.O858 

0.0400 

0.2752 

0.9187 

O.5122 

0.2658 

0.1542 

0.0639 

0.1273 

0.4004 

1.1373 

0.5528 

0.2441 

0.1097 

0.0709 

0.1745 

0.4664 

1.2631 

0.5291 

0.1754 

0.0340 

0.0895 

O.1818 

0.4699 

1.2915 

0.4246 

0.0147 

0.0344 

0.0759 

0.1444 

0.3958 

1.1758 



62 


MCA TN 3014 


TABLE V.- Concluded 

AEBODYNAMIC-INFLUENCE-COEFFICIENT MATRICES 
(b) Concluded 
Plan form 3^1 


0.0954 

-0.0093 

0.0159 

-0.0030 

0.0072 

-0.0010 

0.0023 

0.0102 

0.2938 

0.0514 

0.0266 

0.0133 

0.0103 

0.0040 

0.0022 

0.0564 

0.5381 

0.1427 

0.0526 

0.0317 

0.0135 

0.0020 

0.0138 

0.1356 

0.8003 

0.2488 

0.0854 

0.04l4 

0.0008 

0.0106 

0.0395 

0.2357 

1.0519 

0.3465 

0.1005 

0.0007 

0.0043 

0.0254 

0.0744 

0.3349 

1.2504 

0.3759 

0.0002 

0.0027 

0.0092 

0.0366 

0.0955 

0.3716 

1.2729 



Plan form 3^2 



0.1953 

0.1132 

0.0447 

O.O277 

0.0144 

0.0092 

0.0040 

O.O586 

0.3868 

0.1593 

0.0591 

0.0368 

0.0182 

0.0093 

O.OI67 

0.1136 

0.5389 

0.1934 

0.0688 

0.0404 

0.0157 

0.0084 

0.0344 

0.1589 

O.6818 

0.2311 

0.0776 

0.0365 

0.0037 

0.0190 

0.0496 

0.2059 

0.8162 

0.2663 

0.0744 

0.0023 

0.0084 

0.0274 

0.0642 

0.2502 

0.9230 

0.2655 

0.0009 

0.0043 

0.0100 

0.0296 

O.O679 

0.2594 

0.9282 



Plan form 3^3 



0.1884 

0.1274 

0.0588 

0.0361 

0.0180 

0.0111 

0.0045 

0.0652 

0.3937 

0.1880 

0.0759 

0.0450 

0.0213 

0.0106 

0.0212 

0.1308 

0.5425 

0.2152 

0.0787 

0.0441 

0.0166 

0.0103 

0.0425 

0.1716 

0.6586 

0.2326 

0.0771 

0.0352 

0.0045 

0.0218 

0.0547 

0.2018 

0.7548 

0.2440 

0.0658 

0.0025 

0.0094 

0.0281 

0.0616 

0.2246 

0.8254 

0.2285 

0.0009 

0.0045 

0.0100 

0.0270 

0.0582 

0.2203 

0.8192 



Plan form 344 



0.1799 

0.1332 

0.0699 

0.0435 

0.0217 

0.0125 

0.0048 

0.0678 

0.3917 

0.2080 

0.0917 

0.0528 

0.0240 

0.0112 

0.0245 

0.1431 

0.5450 

0.2352 

0.0893 

0.0473 

0.0169 

0.0119 

0.0497 

0.1842 

0.6459 

0.2374- 

0.0773 

0.0335 

0.0051 

0.0241 

0.0596 

0.2011 

0.7106 

0.2261 

0.0574 

0.0026 

0.0099 

0.0283 

0.0592 

0.2027 

0.7439 

0.1935 

0.0009 

0.0043 

0.0096 

0.0240 

0.0487 

0.1816 

0.7195 
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0.19509 

O.38268 

0.55557 

0.70711 

0.83147 

O.92388 

O.98079 

1.00000 


0.04877 

0.13889 

0.20787 

0.24520 

0.24520 

0.20787 

0.13889 

0.04877 


1.00000 

0.98079 

0.92388 

0.83147 

O.707H 

0.55557 

0.38268 

0.19509 


0.06250 


TABLE VII 

MATRICES USED IN LIFT-DISTRIBUTION CALCULATIONS 


(a) Symmetrical distributions 



0.55557 

0.83147 

0.98079 

0.98079 

0.83147 

0.55557 

0.19509 

0.92388 

0.92388 

0.38268 

-O.38268 

-0.92388 

-0.92388 

-0.38268 

0.98079 

0.19509 

-0.83147 

-0.83147 

0.19509 

0.98079 

0.55557 

0.70711 

-0.70711 

-0.70711 

0.70711 

0.70711 

-0.70711 

-O.7O7H 

0.19509 

-0.98079 

0.55557 

0.55557 

-0.98079 

0.19509 

0.83147 

-0.38268 

-0.38268 

0.92388 

-0.92388 

0.38268 

0.38268 

-0.92388 

-0.83147 

0.55557 

-0.19509 

-0.19509 

0.55557 

-0.83147 

0.98079 

- 1.00000 

1.00000 

- 1.00000 

1.00000 

- 1.00000 

1.00000 

- 1.00000 



[j3in n9 

J - 1 - 




0.09567 

0.13889 

O.I7678 

0.20787 

0.23097 

0.24520 

0.12500" 

0.23097 

0.24520 

O.I7678 

0.04877 

-0.09567 

-0.20787 

-0.12500 

0.23097 

0.04877 

-O.I7678 

-0.24520 

-0.09567 

0.13889 

0.12500 

0.09567 

-0.20787 

-O.I7678 

0.13889 

0.23097 

-0.04877 

-0.12500 

-0.09567 

-0.20787 

O.I7678 

0.13889 

-0.23097 

-0.04877 

0.12500 

-0.23097 

0.04877 

O.I7678 

-0.24520 

0.09567 

0.13889 

-0.12500 

-0.23097 

0.24520 

-O.I7678 

0.04877 

0.09567 

-0.20787 

0.12500 

-0.09567 

0.13889 

-O.I7678 

0.20787 

-0.23097 

0.24520 

-0.12500 



3.00000 5.00000 

2.49441 2.77785 

1.14804 -1.91340 
-O.58527 -4.90395 
-2.12133 -3.53555 
-2.94237 0.97545 

-2.77164 4.61940 

-1.66671 4.15735 


7.00000 

1.36563 

-6.46716 

-3.88899 

4.94977 

5.82029 

-2.67876 

-6.86553 


9.00000 

-1.75581 

-8.31492 

5.00013 

6.36399 

-7.48323 

-3.44412 

8.82711 


11.00000 

-6.11127 

-4.20948 

IO.78869 

-7.77821 

-2.14599 

10.16268 

-9.14617 


13.00000 
-10.80911 
4.97484 
2.53617 
-9.19243 
12 . 75027 
-12.01044 
7.22241 


15.00000 

-14.71185 

13.85820 

-12.47205 

10.60665 

-8.33355 

5.74020 

-2.92635 



0.12500 

0.12500 


0.12500 

0.12500 


0.12500 


0.12500 


0.12500 
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TABLE VII.- Continued 


MATRICES USED IN LIFT-DISTRIBUTION CALCULATIONS 


(a) Concluded 


sln n9 m 

sin 0 m 


1*1 = 


1.00000 

8.54313 

21.30990 

35.19159 

45.24633 

46.88178 

37.02023 

15.00000 

1.00000 

7.24272 

12.07120 

7.00000 

-9.00000 

-26.55640 

-31.38512 

-15.00000 

1.00000 

5.29620 

1.75580 

-10.47641 

-13.46967 

3.86276 

22.95020 

15.00000 

1.00000 

3.00000 

-5.00000 

-7.00000 

9.00000 

11.00000 

-13.00000 

-15.00000 

1.00000 

0.70389 

-5.89795 

4.67719 

6.01353 

-12.97549 

3.05019 

15.00000 

1.00000 

-1.24263 

-2.07105 

7.00000 

-9.00000 

4.55631 

5.38473 

-15.00000 

1.00000 

-2.54328 

2.83220 

-1.39237 

-1.79019 

6.23084 

-11.02088 

15.00000 

1.00000 

-3.00000 

5.00000 

-7.00000 

9.00000 

- 11.00000 

13.00000 

-15.00000 


5.1258 

- 1.8481 

0 

1 

0 

tn 

£ rz 1 

& 

0 

- 0.0481 

0 

-0.0163 

-0.9422 

2.6131 

-1.0193 

0 

-0.1017 

0 

-o.o 4 n 

0 

0 

-0.7022 

1.7999 

-0.7191 

0 

-0.0815 

0 

-0.0226 

- 0.0417 

0 

-0.5649 

1.4142 

-0.5739 

0 

-0.0775 

0 

0 

- 0.0468 

0 

- 0.4881 

1.2027 

-0.4994 

0 

-0.0506 

- 0.0102 

0 

-0.0490 

0 

-0.4494 

1.0824 

- 0.4814 

0 

0 

-0.0160 

0 

-0.0557 

0 

-0.4535 

1.0196 

-0.4106 

-0.0063 

0 

-0.0251 

0 

- 0.0842 

0 

-0.8053 

1.0000 



0.039276 

-0.018861 

- 0.027584 

0 . 01101+8 

- 0.006473 

0.004577 

- 0.003666 

O.OO3236 

-0.00155!+ 


- 0.018861 

0.050969 

- 0.007813 

- 0.034057 

0.015625 

-0.010138 

0.007813 

-O.OO6775 

0.003236 


0.011693 

- 0.026674 

0.044496 

- 0.003236 

- 0.037723 

0.018860 

- 0.013246 

0.011048 

- 0.005220 


-O.OO7813 

0.016913 

- 0.022097 

0.040830 

0 

-0.040830 

0.022097 

-0.016913 

0.007813 


0.005220 

- 0.011048 
0.013246 
-0.018861 
0.037723 
0.003236 
- 0. 044496 
0.026674 
- 0.011693 


-O.OO3236 

O.OO6775 

-O.OO7813 

0.010138 

- 0.015625 

0.034056 

O.OO7813 

- 0.050969 

0.018861 


0.001554 

- 0.003236 

O.OO3666 

- 0.004577 

0.006473 

- 0.011048 

0.027584 

0.018861 

- 0.039276 
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TABLE VII.- Continued 

MATRICES USED IN LIFT-DISTRIBUTION CALCULATIONS 
(b) Antisymmetries! distributions 



0.38268 

0.70711 

0.92388 

1.00000 

0.92588 

0.70711 

0.38268 

0.70711 

1.00000 

0.70711 

0.00000 

-0.70711 

- 1.00000 

-0.70711 

0.92388 

0.70711 

-O.38268 

- 1.00000 

-0.38268 

0.70711 

0.92388 

1.00000 

0.00000 

- 1.00000 

0.00000 

1.00000 

0.00000 

- 1.00000 

0.92388 

-0.70711 

-O.38268 

1.00000 

-0.38268 

-0.70711 

0.92388 

0.70711 

- 1.00000 

0.70711 

0.00000 

-0.70711 

1.00000 

-0.70711 

0.38268 

-0.70711 

0.92388 

- 1.00000 

0.92388 

-0.70711 

0.38268 




[sin n0^] 

= 



0.09567 

0.17678 

0.23097 

0.25000 

0.23097 

0.17678 

0.09567 

0.17678 

0.25000 

O.I7678 

0.00000 

-0.17678 

-0.25000 

-0.17678 

0.23097 

0.17678 

- 0.09567 

-0.25000 

-0.09567 

0.17678 

0.23097 

0.25000 

0.00000 

- 0.25000 

0.00000 

0.25000 

0.00000 

-0.25000 

0.23097 

-0.17678 

- 0.09567 

0.25000 

-0.09567 

-0.17678 

0.23097 

0.17678 

-0.25000 

O.I7678 

0.00000 

-0.17678 

0.25000 

-0.17678 

0.09567 

-0.17678 

0.23097 

-0.25000 

0.23097 

-0.17678 

0.09567 




[cos nej |_n"] 

= 



2.00000 

4.00000 

6.00000 

8.00000 

10.00000 

12.00000 

14.00000 

1.84776 

2.82844 

2.29608 

0.00000 

-3.82680 

-8.48532 

-12.95432 

1 . 41422 

0.00000 

- 4.24266 

- 8.00000 

-7.07110 

0.00000 

9.89954 

O.76536 

- 2.82844 

-5.54348 

0.00000 

9.23880 

8.48532 

-5.35752 

0.00000 

- 4.00000 

0.00000 

8.00000 

0.00000 

- 12.00000 

0.00000 

- 0.76556 

- 2.82844 

5.54328 

0.00000 

-9.23880 

8.48532 

5.35752 

- 1 . 1+1422 

0.00000 

4.24266 

- 8.00000 

7.07110 

0.00000 

-9.89954 

- 1.84776 

2.82844 

-2.29608 

0.00000 

3.82680 

-8.48532 

12.93^32 

- 2.00000 

4.00000 

- 6.00000 

8.00000 - 

10.00000 

12.00000 

- 14.00000 



0.06250 


0.12500 


0.12500 


0.12500 


0.12500 


0.12500 


0.12500 


0.12500 


0.06250 
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TABLE VII. - Concluded 


MATRICES USED IN LIFT- DISTRIBUTION CALCULATIONS 


(b) Concluded 


sin n0„ 


sin 9- 


l?l- 


3.92312 

14.49812 

28.41396 

41.00672 

47.35660 

43.49436 

27.46184 

3.69556 

10.45260 

11.08668 

0.00000 

-18.47780 

-31.35780 

-25.86892 

3.32594 

5.09116 

-4.13292 

-14.39992 

-6.88820 

15.27348 

23.28158 

2.82842 

- 0.00000 

-8.48526 

0.00000 

14.14210 

0.00000 

-19.79894 

2.22228 

-3.40172 

-2.76150 

9.62152 

-4.60250 

-10.20516 

15.55596 

1.53074 

-4.32956 

4.59222 

0.00000 

-7.65370 

12.98868 

-10.71518 

O.78036 

-2.88384 

5.65188 

-8.15672 

9.41980 

-8.65152 

5.46252 


5.1259 

-1.8448 

0.0000 

1 

0 

& G7 1 

** l£j 
11 

0.0000 

-0.0326 

0.0000 

-0.9405 

2.6131 

-1.0138 

0.0000 

-0.0898 

0.0000 

-0.0166 

0.0000 

-O.6983 

1.7999 

-0.7097 

0.0000 

-0.0619 

0.0000 

-0.0396 

0.0000 

-0.5576 

1.4142 

-0.5576 

0.0000 

-0.0396 

0.0000 

-0.04l4 

0.0000 

-0.4742 

1.2027 

-0 . 4666 

0.0000 

-0.0069 

0.0000 

-0.0372 

0.0000 

-0.4199 

1.0824 

-0.3896 

0.0000 

-0.0065 

0.0000 

-0.0286 

0.0000 

-0.3670 

1.0196 


[»]- 


0.040046 

-0.020415 

-0.025984 

0.009366 

-0.004666 

0.002593 

-0.001430 

0.000644 


-0.020415 

0.054108 

-0.011048 

-0.030649 

0.011959 

-0.006097 

O.OO3236 

-0.001430 


0. 014062 
-0.031464 
0.049442 
-0.008456 
-0.032081 
0.012603 
-0.006097 

O.OO2593 


-0.011048 

0.023459 

-0.028872 

0.048012 

-0.007813 

-0.032081 

0.011959 

-0.004666 


0.009396 

-0.019505 

0.022028 

-0.028228 

0.048012 

-0.008456 

-0.030649 

0.009366 


-0.008456 

O.OI7362 

-0.018861 

0.022028 

-0.028871 

0.049442 

-0.011048 

-0.025984 


0.007966 

-0.016269 

0.017362 

-0.019505 

0.023459 

-0.031464 

0.054108 

-0.020415 


-0.003907 

O.OO7966 

-0.008456 

0.009396 

-0.011048 

0.014062 

-0.020415 

0.040046 
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TABLE X 

INTEGRATING MATRICES FOR LOAD AND MOMENT COEFFICIENTS 


L\l = 

|o. 01915 0.03757 0.05454 0.06942 0.08163 0.09070 0.09629 0.04909J 

[^ c bm| ' 

[0.01878 0.03472 0.04533 0.04913 0.04528 0.03487 0.01833 0.00196J 

^ Cli/2 J 

[0.01993 0.03596 0.05713 O.06556 0.08736 0.08164 O.II365J 




00939 0.01735 


N ‘ 

0.02267 0.02454 0.02267 0.01736 O.OO939J 
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MCA TN 3014 



(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 1.- Spanwise lift distributions for wings of very low aspect 

ratio ( from ref . 3 ) • 


Spanwise loading coefficient, — Spanwise loading coefficient, =g 
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Dimensionless spanwise ordinate, y* 


(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 
Figure 1.- Concluded. 
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Dimensionless spanwise ordinate, y* 


(a) Symmetrical lift distributions. 

(b ) Anti symmetrical lift distributions. 

Figure 2.- Spanwise lift distributions for plan form 3 11 (A = 

A = 0). 



loading coefficient, — Spanvise loading coefficient 
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0 •! .2 .3 .4 .5 .6 .7 .8 .9 1.0 


Dimensionless spanwlse ordinate, y* 

(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 


Figure 2.- Concluded. 
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(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 3.- Spanwise lift distributions for plan form 312 (A = 1.5j 

A = 0.25). 


Spanwise loading coefficient, — Spanwise loading coefficient, 
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Dimensionless spanwise ordinate, y* 

(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 


Figure 3*- Concluded. 
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(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 4.- Spanwise lift distributions for plan form 313 (A = 1.5} 

A = 0.50). 


NACA TN 5014 


87 



0 <1 .2 .3 .4 .6 .6 .7 .8 .9 1.0 


Dimenalonlesa apanwlae ordinate, y* 

(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 
Figure 4.- Concluded. 
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Dimensionless spanwise ordinate, y* 


(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 5 .- Spanwise lift distributions for plan form 3 l 1 

A = 1.00). 


Spanwise loading coefficient, =p — Spanwise loading coefficient 
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0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

Dimensionless spanwise ordinate, y* 


(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 


Figure 5*- Concluded. 




Spanwise loading coefficient, Spanwise loading coefficient 
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(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

- Spanwise lift distributions for plan form 315 (A = 1.5; 

A = 1.50). 


Figure 6. 



Spanwise loading coefficient, =p — Spanwise loading coefficient 
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Dimensionless spanwise ordinate, y* 

(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 
Figure 6.- Concluded. 
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(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure r J.~ Spanwise lift distributions for plan form 321 (A = 3 0* 

A = 0). ’ 



Spanwise loading coefficient, ^ — S panwise loading coefficient. 
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Dimensionless spanwise ordinate , y* 


(c) Lift distribution for inboard flap. 

(d ) Lift distribution for outboard aileron. 


Figure 7*- Concluded. 
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(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 8.- Spanwise lift distributions for plan form 322 (A = 3.0; 

A = 0.25). 




Spanwise loading coefficient, -=g — Spanwise loading coefficient 
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Dimensionless spanwise ordinate, y* 


(c) Lift distribution for inboard flap. 

(d ) Lift distribution for outboard aileron. 

\ 

Figure 8.- Concluded. 
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(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 9.- Spanwise lift distributions for plan form 323 (A = 3.0; 

A = 0.50). 



Spanwise loading coefficient, — Spanwise loading coefficient, 
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Dimensionless spanwise ordinate, y* 


(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 


Figure 9« - Concluded. 






Spanwise loading coefficient, Spanwise loading coefficient 
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(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 


Figure 10.- Concluded. 



Spanwise loading coefficient, Spanwise loading coefficient, - 
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Dimensionless spanwise ordinate, y* 


(a) Symmetrical lift distributions. 

(b ) Antisymmetrical lift distributions. 

Figure 11.- Spanwise lift distributions for plan form 325 (A = 3-0j 

A = 1.50). 



Spanwise loading coefficient, — Spanwise loading, coefficient, 
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Dimensionless spanwise ordinate, y* 


(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 


Figure 11.- Concluded. 
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Dimensionless spanwise ordinate, y* 


(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 12.- Spanwise lift distributions for plan form 331 (A = 6 . 0 ; 

A = 0). 



Spanwise loading coefficient, — i- Spanwise loadlng coelflcient) _ f 
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(c) Lift distribution for inboard flap. 

(d ) Lift distribution for outboard aileron. 


Figure 12.- Concluded. 
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(a) Symmetrical lift distributions. 


(b) Antisymmetrical lift distributions. 

Figure 13.- Spanwise lift distributions for plan form 332 (A = 6.0; 

A = 0.25). 


Spanwise loading coefficient, =£ — Spanwise loading coefficient, =£ 
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(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 
Figure 13 . - Concluded. 
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(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 14.- Spanwise lift distributions for plan form 333 (A = 6.0; 

A = 0.50). 



Spanwise loading coefficient, Spanwise loading coefficient, 
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Dimensionless spanwise ordinate, y* 


(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 
Figure 14.- Concluded. 
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Dimensionless spanwise ordinate, y* 


(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 15 .- Spanwise lift distributions for plan form 334 (A = 6.0; 

A = 1.00). 





Spanwise loading coefficient, — Spanwise loading coefficient 
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Dimensionless spanwise ordinate , y+ 


(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 


Figure 15 . - Concluded. 
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(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 1 6 .- Spanwise lift distributions for plan form 355 (A = 6.0; 

A = 1.50). 


Spanwise loading coefficient, — Spanwise loading coefficient, 
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Dimensionless spanwise ordinate, y* 

(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 
Figure l6. - Concluded. 
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(a) Symmetrical lift distributions. 

(b) Antisymmetrical lift distributions. 

Figure 17. - Spanwise lift distributions for plan form 3*+l (A = 12.0; 

A = 0). 


Spanwise loading coefficient, - — — Spanwise loading coefficient, r 

cC, , c 
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Dimensionless spanwise ordinate, y+ 


( c ) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 
Figure 17* - Concluded. 
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Figure 18. - Spa 
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Dimensionless spanwise ordinate, y* 

(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 
Figure 18 . - Concluded. 
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Dimensionless s panwise ordinate, y* 


(c) Lift distribution for inboard flap, 
(d) Lift distribution for outboard aileron. 


Figure 19. - Concluded. 
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Dimensionless spanwise ordinate, y* 


(a) Sy m metrical lift distributions. 

(b) Anti symmetrical lift distributions. 

Figure 20.- Spanwise lift distributions for plan form 344 (A = 12. 0; 

A = 1.00). 
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Effective angle of attack for unit flap deflection, 
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Flap - chord ratio, Cffc 

Figure 21.- Variation of flap effectiveness with flap-chord ratio. 


NACA- Langley - 9-30-53 - 1000 


